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Summary  Narrative. 


In  this  program  the  measurement  of  stresses  in  ceramics  using  their  optical 
fluorescence  is  being  developed.  The  technique  is  based  on  the  phenomena  of  piezo¬ 
spectroscopy,  namely  the  shift  in  optical  fluorescence  of  a  material  produced  by  a 
stress.  In  the  course  of  the  program  to  date  we  have  a)  derived  for  the  first  time  a 
general  relationship  between  an  observed  shift  in  fluorescence  frequency  and  the 
state  of  stress,  b)  applied  it  to  the  determination  of  both  applied  and  residual  stresses 
in  single  crystal  and  polycrystalline  materials,  and  c)  begun  to  exploit  the  techruque 
to  the  measurement  of  stresses  in  a  variety  of  composite  materials.  These  are 
described  in  detail  in  the  accompanying  technical  reports  appended  here. 

As  part  of  this  investigation  we  have  also  found  that  in  ruby  (Cr^"^  doped  sapphire) 
the  observed  optical  fluorescence  shift  as  a  function  of  chromium  concentration  is  a 
direct  measure  of  the  "image  stress".  This  component  of  stress,  originating  from 
regions  of  mis-fit  with  the  matrix,  was  introduced  by  Eshelby  in  the  early  1950’s  but 
has  not  hitherto  been  measured  directly.  Although  no  application  is  yet  evident  to 
us,  the  measurement  of  the  "image  stress"  does  have  an  important  place  in  the 
conceptual  framework  of  the  elasticity  of  inclusion  and  composite  theory. 

On-going  work  that  has  not  yet  reached  completion  includes: 

a)  development  of  a  process  for  incorporating  two  different  dopant  species,  each  of 
which  fluoresces  at  a  different  frequency.  This  is  desirable  since  simultaneous 
fluorescence  at  different  frequencies  would  facilitate  the  determination  of  the 
complete  stress  tensor.  We  have  successfully  demonstrated  simultaneous 
fluorescence  from  two  co-dopants  but  need  now  to  measure  the  piezo-spectroscopic 
coefficients  for  the  different  dopants. 

b)  determination  of  the  residual  stresses  in  two-phase  composites  as  a  function  of 
the  volume  fraction.  Using  the  methodology  developed  in  the  first  year,  and 
described  in  report  #1  attached,  we  have  made  measurements  of  the  average 
residual  stresses  and  the  second  moment  of  the  stress  distribution  for  a  range  of 
AlaOa-ZrOa  composites.  We  are  now  using  an  "information  theory"  approach  to 


relate  these  measurements  to  the  thermal  expansion  coefficients  and  derive  the 
temperature  at  which  the  thermal  expansion  stresses  began  to  develop  on  cooling. 

c)  measurement  of  the  stress  distribution  along  embedded  fibers  in  a  composite 
from  the  shape  of  the  fluorescence  peak  generated  by  a  light  beam  propagating  along 
and  through  an  embedded  fiber  in  a  composite. 
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Abstract 

A  general  methodology  is  developed  for  determining  the  state  of  stress  and  the 
numerical  value  of  the  stresses  from  observed  shifts  and  broadening  of  optical 
fluorescence  lines.  The  method  is  based  on  the  piezo-spectroscopic  properties  of 
single  crystals.  We  present  general  relationships  between  the  measured  fluorescence 
shifts  and  the  stress  state  for  a  number  of  illustrative  cases,  pertinent  to  both 
polycrystalUne  and  single  crystal  ceramics  under  stress.  These  include  measuring  the 
stresses  applied  to  poy crystalline  ceramics,  the  residual  stress  distribution  due  to 
crystallographically  anisotropic  thermal  expansion,  and  the  stresses  applied  to  single 
crystals.  Using  the  recently  implemented  technique  of  performing  the  fluorescence 
measurements  in  an  optical  microprobe,  we  also  provide  experimental  tests  of  the 
relationships  derived. 


1  INTRODUCTION 


Despite  the  importance  of  measuring  the  stress  state  in  a  material  or  structure  there  are 
few  quantitative  methods  available.  Those  that  have  been  developed  and  widely  used,  such  as 
photoelastic  methods  and  X-ray  methods,  whilst  accurate  are  generally  limited  to  examining 
flat  surfaces  and  ones  where  the  area  interrogated  is  not  obscured  or  shadowed  by  features 
such  as  steps.  In  addition  they  are  difficult  to  apply  to  the  measurement  of  stresses  from 
small  areas,  a  few  microns  across,  and  are  not  generally  microscope- based  thus  precluding 
making  measurements  from  features  identifiable  during  examination  of  a  microstructure. 

For  these,  and  related,  reasons  a  number  of  optical  techniques  based  on  piezo-spectroscopic 
phenomena  have  been  proposed.  One  gaining  increased  support,  especially  in  VLSI  silicon 
based  technologies,  is  that  of  Raman  spectroscopy  using  an  optical  microprobe  [1].  The 
Raman  lines  of  silicon  at  520  cm*'  are  shifted  by  the  application  of  a  stress,  and  so  re¬ 
gions  on  a  silicon  wafer  around  fine  metal  lines  and  dielectric  features  can  be  probed,  with 
approximately  a  micron  resolution,  and  their  stress  evaluated,  .\nother  piezo-spectroscopic 
technique  proposed  in  the  late  lOTO’s  for  the  measurement  of  stress,  is  that  of  fluorescence 
spectroscopy  [2|.  The  suggestion,  made  by  Grabner,  followed  the  demonstration  at  the  Bu¬ 
reau  of  Standards  that  the  pressure  in  high-pressure  experimental  cells  could  be  calibrated 
by  the  shift  of  the  chromium  fluorescence  lines  from  fine  particles  of  ruby  embedded  in  the 
high-pressure  cell  [3j.  In  turn,  this  work  built  on  the  luminescence  studies  of  chromium 
doped  sapphire  by  Schawlow  [4]  and,  independently  by  Sugano  and  Tanabe  [5],  as  part  of 
their  investigation  of  possible  laser  materials.  They  found  that  the  luminescence  lines  in 
chromium-doped  sapphire  were  shifted  by  an  applied  stress.  VVe  have  recently  demonstrated 
that  by  implementing  the  fluorescence  spectroscopy  using  an  opticaj  microprobe  it  is  pos- 
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sible  to  measure  stresses  in  chromium-doped  alumina  from  regions  identified  in  the  optical 
microscope  and  with  a  spatial  resolution  of  a  few  microns  [6]. 

In  this  paper  we  present  a  method  of  analyzing  an  observed  frequency  shift  in  terms  of 
the  piezo-spectroscopic  coefficients  and  the  stress  (and  hence  strain)  on  the  crystal  structure 
containing  a  fluorescing  atom.  This  provides  a  correction  to  the  equation  derived  by  Grab- 
ner  for  the  case  of  chromium  fluorescence  in  alumina  [2).  VV'e  then  extend  the  analysis  to 
the  measurement  of  stresses  in  inhomogeneously  stressed  materials,  presenting  a  number  of 
illustrative  examples  of  the  application  of  the  fluorescence  technique  to  the  measurement  of 
commonly  occurring  types  of  stresses  in  ceramics.  These  provide  the  basis  for  future  work  in 
the  widespread  application  of  the  technique  to  measure  stresses  in  ceramic  composites  and 
at  ceramic  interfaces. 

2  PRINCIPLES  OF  PIEZO-SPECTROSCOPIC  ANALYSIS 

The  basis  of  all  piezo-spectroscopic  methods  for  the  measurement  of  stress  in  crystalline 
materials  is  that  an  applied  stress  strains  the  lattice  and  alters  the  energy  of  transitions 
between  electronic  or  vibrational  states.  In  this  way.  the  energy  of  any  radiative  transitions, 
such  as  those  that  give  rise  to  luminescence,  also  varies  and  produces  systematic  shifts 
in  the  frequency  of  luminescence  lines  such  as  fluorescence  lines.  In  common  with  other 
spectroscopies,  the  sharpness  of  the  fluorescence  line  depends  on  temperature  and  variations 
in  the  site  position  of  the  fluorescing  atomic  species.  When  the  fluorescing  ion  is  located  on 
a  lattice  site  in  a  tightly  bound,  high  elastic  modulus  material,  eis  is  the  case  for  chromium 
substitution  for  an  aluminum  ion  in  the  sapphire  lattice,  the  lines  can  be  exceptionally  sharp, 
broadened  only  by  phonon  scattering.  Interstitial  dopants,  where  there  is  greater  variability 
in  the  dopant  site  or  in  a  random  network  structure,  such  as  chromium  in  glaisses,  give  rise 
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to  much  broader,  and  usually  weaker,  luminescence  lines. 

The  principle  of  relating  an  observed  line  shift  in  a  fluorescence  or  absorption  spectrum 
to  the  state  of  stress  have  been  described  previously  by  Grabner  [2J  in  reference  to  the 
measurement  of  stresses  in  sintered,  polycrystalline  alumina  containing  a  trace  of  Cr^"^.  The 
more  general  case  is  presented  here. 

When  a  crystal  is  subject  to  a  homogeneous  stress  a,’  (in  crystallographic  frame  of  refer¬ 
ence),  the  change  in  frequency.  Af/,  of  a  luminescence  line  is  given  by  the  tensorial  relation¬ 
ship: 

Ai/  =  n„<T,*  (I) 

where  H.^'s  are  the  piezo-spectroscopic  coefficients  relating  frequency  to  stress.  They  form 
symmetrical  second  rank  tensors,  which  have  symmetry  properties  governed  by  the  point 
symmetry  group  of  the  fluorescing  ion.  They  may  also  have  different  numerical  values  for 
different  fluorescence  line  energies.  No  appreciable  stress  induced  broadening  has  been  re¬ 
ported  at  room  temperature  so  the  lines  merely  shift  with  stress  without  changing  shape. 
(.Stress  induced  splitting  is  detectable  at  liquid  helium  temperatures  but  is  unobservable  at 
room  temperature  being  masked  by  thermal  broadening).  In  general,  any  applied  stress  will 
not  be  imposed  on  the  crystal  along  its  crystallographic  a.xes,  as  represented  in  equation  1, 
but  rather  in  some  other  coordinate  frame.  The  stress  on  the  crystal  structure  will  then  be 
given  by  the  applied  stress  components  resolved  onto  the  crystallographic  a.xes,  i.e.  through 
an  orthogonal  coordinate  transformation: 

=  OifcUj/ (2) 

where  is  the  transformation  matrix  relating  the  frame  of  reference  in  which  the  applied 
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stresses  are  defined  with  respect  to  the  crystallographic  axes  of  the  lattice. 

Thus,  the  frequency  shift  of  a  fluorescence  line  in  a  luminescing  crystal  oriented  at  an 
arbitrary  angle  to  a  superimposed  stress  and  strain  field  is  given  by  the  follovung  tensorial 
relation; 


Ai/  =  riij  aiitOjl  (Tkl  =  n,j  a/cmflln  Cn 


(3) 


One  of  the  major  limitations  of  piezo- spectroscopy  is  evident  from  this  equation.  Unless 
the  fluorescing  species  emits  at  several  different  frequencies,  so  that  a  set  of  simultaneous 
equations  (one  for  each  transition)  can  be  generated,  there  may  be  an  insufficient  number  of 
distinct  transitions  to  solve  for  each  of  the  maximum  of  six  independent  stress  components 
needed  to  completely  describe  the  stress  state.  In  some  materials,  the  point  symmetry  of  the 
fluorescing  ion  may  also  preclude  the  determination  of  all  the  independent  stress  components. 
This  will  be  seen  to  be  the  case  for  both  MgO  and  .\l2O3. 

Of  the  many  fluorescence  spectra  reported,  that  of  chromium  fluorescence  in  sapphire 
has  probably  been  the  most  thoroughly  studied.  The  chromium  ions  replace  aluminum  ions 
substitutionally  in  the  sapphire  lattice  and  maintain  their  trivalent  state,  Cr^'*'.  .Although 
the  chromium  ion  is  surrounded  by  an  octahedron  of  oxygen  ions,  the  octahedron  is  distorted 
along  the  body  diagonal  giving  a  point  symmetry  of  D.^dlJm)  for  the  chromium  ion.  The 
distortion  also  causes  an  asymmetry  of  the  electrostatic  crystal-field  and  results  in  a  splitting 
of  the  energy  levels  by  0.004  eV,  in  turn  giving  rise  to  the  two  distinct  radiative  transition 
lines  in  ruby,  Ri  and  Rj  at  1.790  and  1.794  eV  respectively.  From  the  point  symmetry  of  the 
chromium  ion,  it  can  be  shown  that  the  off-diagonal  components  of  the  piezo-spectroscopic 
tensor  are  zero  and  that  the  diagonal  components  are  symmetrical  about  the  c-axis,  namely 
fill  =  U.22  ^  1133.  For  simplicity,  we  will  use  IIu  =  1122  =  Oa  and  Ilaa  =  Re  for  Ri  and  R2 
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lines.  Substituting  these  components  into  equation  3,  the  frequency  shift  for  a  chromium 
fluorescence  line  can  be  expressed  in  terms  of  the  applied  stress  components  as: 

Ai/  =  nalO"!!  4-  (722  +  <733)  +  (Oc  —  na)(a3j  <7ii  032(722  +  033(733) 

-i-  2(nc  —  na)(a3ia32  <712  +  O32O33  (T23  +  031033  <731)  (4) 

This  differs  from  the  expression  derived  by  Grabner  (equation  4  reference  2)  in  having  an 
additional  term  (the  second)  proportional  to  the  trace  of  the  applied  stress  tensor.  For  the 
particular  case  in  which  the  applied  stress  is  purely  hydrostatic,  this  equation  simplifies  to 
that  used  in  the  high-pressure  physics  and  materials  literature,  namely 

Ai/ = -{2na  +  n,)P  (5) 

where  P  is  the  pressure.  .Also,  unlike  Grabner’s  equation,  equation  4  reduces  to  equation  1 
when  the  stress  field  is  applied  along  the  crystal  axes.  The  point  symmetry  of  the  chromium 
ion  simplifies  equation  3  but  nevertheless  the  stress  components  cannot  be  determined  com¬ 
pletely  from  the  shifts  of  the  Ri  and  R2  lines  unless  the  misorientation  between  the  applied 
stress  and  the  crystal  axes  is  known  from  some  other,  independent  measurement. 

The  chromium  fluorescence  spectrum  from  magnesia  heis  also  been  investigated  in  some 
detail  [7],  and  can  similarly  be  used  to  monitor  strains.  As  in  sapphire,  the  chromium  ion 
occupies  substitutionally  a  lattice  site  in  magnesia  and  is  surrounded  by  an  octahedron  of 
oxygen  ions.  However,  the  oxygen  octahedra  remain  undistorted  so  the  chromium  point 
symmetry  is  raised  to  m3m.  As  a  consequence,  the  level  is  not  split  and  only  one  line, 
at  14319  cm"*,  occurs.  The  m3m  point  symmetry  imposes  the  additional  requirement  that 
the  diagonal  piezo-spectroscopic  coefficients  are  all  equal  (Hu  =  022  =  ^33).  Thus,  for 
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chromium  luminescence  in  magnesium  oxide,  the  line  shift  is  sensitive  only  to  the  trace  of 
the  stress  tensor: 

Aj/ =  nn(<7u  +  <^22  +  (^ss)  (6) 

This  insensitivity  to  the  deviatoric  stresses,  nevertheless,  has  the  advantage  that  the  shift  is 
independent  of  the  orientation  of  the  crystal  structure  with  respect  to  the  applied  stress. 

Rare  earth  ions,  such  as  Sm^"*"  and  Nd^"*",  can  substitute  for  the  yttrium  ion  in  the  yt¬ 
trium  aluminum  oxide,  garnet  structure  and  give  rise  to  sharp  cheiracteristic  fluorescence 
spectra.  The  six  possible  substitutional  sites,  although  equivalent  under  the  cubic  point 
group  symmetry  of  YAG,  have  orthorhombic  site  symmetry.  Since  ions  at  all  three  sites  will 
fluoresce  when  the  crystal  is  excited,  applying  a  stress  should  in  general  break  the  degen¬ 
eracy  associated  with  the  equivalence  of  the  three  sites  and  cause  line  splitting.  However, 
in  practice,  the  stress  splitting  is  masked  by  thermal  broadening  at  room  temperature,  and 
so  the  fluorescence  appears  to  come  from  a  crystal  having  cubic  symmetry  [8]  with  equal 
piezo-spectroscopic  coefficients. 

3  INHOMOGENEOUS  STRAIN  FIELDS 

The  above  analysis  pertains  to  the  case  of  a  single  crystal  subject  to  a  homogeneous  strain 
field  over  the  volume  of  material  exposed  to  the  excitation  probe.  When  the  probed  volume 
contains  more  than  one  crystallite  or  the  strains  vary  within  it,  the  observed  spectral  shifts 
are  more  complex.  They  can,  however,  be  understood  from  simple  superposition  arguments 
provided  inel^lstic  scattering  of  the  luminescence  is  not  appreciable  and  the  strains  do  not  vary 
over  distances  comparable  to  the  lattice  spacing.  (For  the  case  of  the  N-  and  F-  fluorescence 
lines,  the  latter  requirement  is  a  little  more  stringent.  In  that  case,  the  strains  should  not 
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vary  over  distances  associated  with  the  separation  of  the  ions  making  up  the  chromium  pairs 
from  which  these  characteristic  lines  originate). 

Subject  to  these  restrictions,  strain  inhomogeneities  and  the  presence  oi  multiple  grams 
within  the  probed  volume  can  cause  both  broadening  and  shifting  of  the  luminescence  lines 
as  analyzed  in  .Appendix  .A.  In  essence,  each  element  of  material  over  which  there  is  a 
constant  crystallographic  orientation  and  homogeneous  strain  contributes  a  signal  intensity 
proportional  to  it. 5  volume  and  concentration  of  fluorescence  ions.  The  spectral  shift  is  given 
by  equation  3  with  the  appropriate  values  for  each  of  the  parameters.  The  superposition  of 
the  lines,  shifted  according  to  the  strain  of  each  volume  element,  causes  the  observed  line 
to  appear  to  be  broadened  and  shifted.  The  broadening  is  not  a  measure  of  the  deviatoric 
strains  per-se  as  has  been  staled  in  the  literature  previously  [2]  but  rather  of  the  statistical 
variation  in  strains  and  crystallographic  orientations  in  the  region  examined.  As  a  result  the 
broadening  need  not  be  symmetrical.  Indeed  the  b-oadened  lines  can  have  a  fine  structure 
resulting  from  there  being  one  or  two  regions  whose  contributions  are  disproportionately 
large.  This  is  probably  the  cause  of  the  fine  structure  in  the  spectra  reported  by  Grabner 

.Although  analytical  relationships  between  line  shifts,  broadening  and  inhomogeneous 
stress  state  in  polycrystalline  materials  can  not,  in  general,  be  presented,  a  number  of  simple 
cases  of  practical  importance  can  be  considered.  Equations  for  a  number  of  these  are  derived 
as  described  in  the  following  sub-sections. 

3,1  Polycrystalline  Material 

The  first  example  for  which  an  analytical  function  can  be  derived  corresponds  to  a  poly¬ 
crystalline  material,  with  no  crystallographic  texture  or  internal  residual  stress,  under  a 
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homogeneous  applied  stress  and  probed  with  an  excitation  beam  large  compared  to  the 
grain  size.  In  this  C'lse,  the  contributions  from  the  individual  crystallites  are  averaged  over 
many,  randomly  oriented  grains. 

The  frequency  shift  in  the  fluorescence  line  is  given  by  the  first  moment  of  the  distribution 
of  fluorescence  frequencies.  Using  equations  3,  and  .A6  and  B3  in  the  appendices,  this  can 
be  simply  expressed  as: 

Ai/  —  Jjj  V')At'dd(i0d0 

which  reduces  to 

J  (^11  +  n22  +  n33)(<Tii  +  0-22  +  <T33)  (7) 

This  equation  shows  that,  irrespective  of  the  crystal  structure,  the  frequency  shift  of  the 
fluorescence  from  a  large  number  of  grains  in  a  poiycrystalliue  material  is  dependent  only 
on  the  trace  of  the  applied  stress  tensor. 

The  polycrystalline  ensemble  of  grains  also  leads  to  an  intrinsic,  stress  dependent  broad¬ 
ening  of  the  fluorescence  line.  The  peak  broadening  is  given  by  the  second  moment  of  the 
distribution,  which  using  equation  A8,  may  be  written  as: 


<  Ai/^  >  =  —  J(n22  —  +  (1133  —  —  ((1122  —  nii)tn33  —  rin)] 

X  ~  <^U<^22  “  (^22^33  “  C’33<7u)  -f  -t-  -f 
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(8) 


Such  broadening  is  due  to  the  differences  in  the  diagonal  piezo-spectroscopic  coefficients. 
Consequently,  for  materials  with  fluorescing  iorts  having  cubic  point  symmetry,  such  as  Cr^"*" 
in  MgO,  (for  which  the  piezo-spectroscopic  coefficients  are  equal)  the  application  of  an  ex¬ 
ternal  stress  does  not  cause  any  broadening. 
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3.2  Inhomogeneously  Stressed  Single  Crystal 

The  second  case  corresponds  to  a  single  crystal  containing  inhomogeneously  stressed  re¬ 
gions.  A  special  case,  but  nevertheless  an  important  one,  is  one  in  which  the  stress  varies 
with  depth  below  the  surface,  such  as  might  pertain  to  surface  stressing  produced  by  ion 
exchange  or  when  the  sample  is  loaded  in  benJ>ng.  For  purposes  of  generality,  we  may  as¬ 
sume  that  the  concentration  of  the  fluorescing  ion,  c,  is  a  function  of  depth  z.  We  can,  in 
principle,  also  incorporate  various  instrumental  functions,  such  as  the  collection  efficiency 
as  a  function  of  depth  of  the  apparatus  used  to  meaisure  the  fluorescence  spectrum,  into  a 
function  g{z).  The  function  g{z)  allows  the  fluorescence  signal  to  be  normalized: 

j  ciz)g{z)dz  =  1 

and  a  weighting  function  introduced: 

W{Su)diy  =  c{z)g{z)dz  (9) 

Again,  for  purposes  of  generality,  the  frequency  shift  ca  be  written  so  as  to  include  a  shift 
dependent  on  the  concentration,  c,  of  the  fluorescing  species; 

=  U.,j  a,kaji  (7ki  +  3{c  -  Co)  (10) 

where  cq  is  the  concentration  of  the  reference  sample  which  the  shift  is  compared. 

The  average  peak  shift,  expressed  by  the  first  moment  of  the  fluorescence  peak  then 
becomes 

Ay  =  j[U,,aikajicrk!  +  co)]c{z)g{z)dz  (11) 
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The  broadening  is  a  mixture  of  the  inhomogeneous  broadening  due  to  the  depth  depen¬ 
dence  of  shift  and  the  intrinsic  homogeneous  broadening  due  to  change  of  concentration.  For 
simplicity,  we  give  the  result  for  the  case  of  constant  concentration  cq: 

<  >=  j{U,ja,kaji(Ty)'^Ci)g{z)dz  -  (12) 

in  which  Su  ^  is  also  evaluated  at  c  =  co- 

3.3  Anisotropic  Thermal  Expansion  Residual  Stresses 

The  third  case  considered  here  is  of  a  polycrystalline  materia]  for  which  the  thermal  ex¬ 
pansion  coefficient  is  crystallographically  anisotropic.  .As  a  result  of  such  thermad  expansion 
anisotropy,  residual  stresses  can  be  created  on  cooling  from  the  fabrication  temperature,  with 
the  magnitude  of  the  stresses  being  dependent  on  the  geometrical  constraint  imposed  by  the 
neighboring  grains,  the  elastic  properties  of  the  grains  and  the  extent  to  which  accommo¬ 
dation  mechanisms,  such  as  creep  and  plastic  deformation,  which  can  relax  the  stresses.  In 
the  special  cases  in  which  the  grains  do  not  constrain  one  another  on  cooling  or  the  accom¬ 
modation  mechanisms  are  complete,  no  internal  stresses  de\.elop  and  neither  a  line  shift  nor 
broadening  would  arise.  However,  there  is  in  general  some  mutual  constraint  imposed  by  the 
contracting  grains  that  causes  stresses  to  develop  within  the  material.  The  average  stress 
over  the  material  must  be  zero  but  variations  in  stress  from  one  grain  to  another  can  cause 
both  a  broadening  of  the  line  and  a  line  shift.  .Although  the  specific  values  will  vary  from 
place  to  place  in  the  sample,  since  the  stresses  will  also  vary  from  grain  to  grain,  it  is  possible 
to  calculate  the  average  values  as  will  now  be  shown.  Two  assumptions  are  made.  First, 
that  the  probed  volume  encompasses  a  sufficiently  large  number  of  grains  that  an  ensemble 
average  can  be  considered,  and  secondly  that  the  distribution  of  stresses  can  be  described 
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by  an  appropriate  distribution  function. 


Consider  a  material,  such  as  alumina,  which  has  a  thermal  expansion  coefficient  that  is 
larger  along  its  c-axis  than  along  its  a-axis,  and  has  an  axis  of  rotational  symmetry  along 
the  c-axis.  On  cooling  from  its  fabrication  temperature,  tensile  stresses  will,  on  average, 
develop  along  the  c-axis  of  the  grains  and  compressional  stresses  will  develop  in  the  basal 
plane.  Therefore, 

5r<0,  ^>0  (13) 

Now,  consider  a  grain  intersected  by  an  arbitrarily  placed  plane  that  divides  the  poly¬ 
crystalline  sample  into  two.  The  orientation  of  the  grain  can  be  specified  by  the  angular 
coordinates  (ff,  0)  of  its  c-axis  (figure  1).  Provided  the  dividing  plane  intersects  a  large  num¬ 
ber  of  grains,  with  all  possible  stress  states,  we  need  only  consider  the  average  stresses  along 
the  a-  and  c-  axes.  Therefore  the  stress  in  the  basal  plane  is  equi-biaxial,  and  rotation  of 
the  grain  around  the  c-axis  does  not  alter  its  stress  state.  The  force  acting  per  unit  area  in 
any  grain  perpendicular  to  the  dividing  plane  can  be  expressed  by  resolving  the  stresses  in 
that  grain: 

/  =  9^sin*5 -t-^cos^^  (14) 


Since,  across  the  dividing  plane,  the  orientation  of  the  grains  is  random,  the  average  force 


can  be  written  as: 


_  1  f2!r  ,f/2 

^  =  27/0  ‘‘*1 


sin 


Since  the  force  must  on  average  be  zero,  /  =  0,  the  mean  stresses  in  the  c-  and  a-  crystallo¬ 
graphic  directions  are  related  by: 

5c  =  -257  (16) 
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Recognize  that  the  average  stresses  described  above  are  equivalent  to  volume  average,  the 


mean  shift  of  the  fluorescence  line  is  therefore, 

A//  =  {lie  - 


(17) 


The  peak  broadening  can  be  calculated  using  equation  A8  on  the  assumption  that  the 
distribution  in  stresses  in  the  a-  and  c-  crystallographic  directions  are  independent  of  one 
another; 

Pi<^a,<Tc)  =  pi<Ta)p{<yc)  (18) 

and  are  also  Gaussian  with  the  same  width,  <Tw 


p{<7a)  —  —  exp 

V2k<t^ 


(19) 


The  assumption  of  independence  of  the  stress  distributions  in  the  different  directions  enables 
the  weighting  function  in  equation  A8  to  be  written  as; 


W(Ai/)di/  =  P((Taj<7c)d<r,dcrc  =  p(crj)p{<Tc)dcrad(7c 


Thus,  using  measurements  of  the  mean  peak  shift,  which  from  equation  17  gives  the  mean 
stresses  along  the  a-  and  c-  axes,  and  a  measurement  of  the  peak  width,  the  stress  distribution 
within  the  material  can  be  determined. 
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These  equations  shows  that  the  crystallographic  anisotropy  in  thermal  expansion  leads 
to  both  an  average  shift  and  broadening  of  the  fluorescence  line.  They  are  related  to  the 
average  values  of  the  stresses  and  the  stress  distribution.  The  spatially  varying  stresses  can 
be  significantly  larger  than  the  average,  especially  those  acting  over  smaller  distances  than 
the  grain  size,  related  to  the  presence  of  the  grain  edges  and  corners.  These  are  primarily 
responsible,  for  instance,  for  micro-crack  extension.  Although  detailed  information  of  such 
stresses  does  not  appear  to  be  amenable  to  the  analysis  of  the  shape  of  the  fluorescence  lines, 
their  magnitude  and  the  occurring  frequency  are  reflected  by  the  stress  distribution. 

3.4  Matrix  Stresses  Due  to  Presence  of  Fibrous  Inclusions 

The  final  case  considered  is  that  of  a  polycrystalline  matrix  surrounding  an  embedded 
fiber.  The  radial  and  tangential  stresses  in  a  matrix  around  a  misfitting  fiber  of  radius,  Tq, 
are  given  by  standard  elastic  solutions,  varying  with  the  inverse  square  of  distance,  r,  from 
the  fiber  axis.  For  observations  made  in  a  cross-section  perpendicular  to  the  fiber  axis  (axis 
3),  the  matrix  stresses  may  be  written  as 


(7ii  =  — <T22 


(Tori  , 

’ 


<733  =  0 


(21) 


Thus,  from  equation  3,  the  frequency  shifts  are  given  by 


At'  =  (1122  ~  ~  <^22)^11  "h  (^33  ~  nil)(<^3l  ~  *^32)  *^22 


which  would  be  detectable  with  a  probe  having  high  optical  resolution  (<C  ro).  If  the  probe 
is  large,  such  that  either  many  grains  in  a  polycrystalline  matrix  or  regions  within  a  circle 
encompassing  the  fiber  in  a  single  crystal  matrix  are  excited  simultaneously,  then  the  average 
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peak  shift  is  At/  =  0.  The  peak  broadening,  for  such  a  large  probe,  would  also  average  out 
for  single  crystal  and  polycrystalline  material  to  a  value  of 

<  >=  ~  (Til  [(022  -  riu)^  +  (033  -  flu)^  -  (022  “  nn)(n33  -  On)]  (23) 

where  is  averaged  over  the  volume  probed. 

4  OBSERVATIONAL  TESTS 

In  this  section  a  series  of  measurements  are  described  that  provide  experimental  tests  and 
validations  of  the  equations  presented  in  the  preceding  sections.  The  tests  were  made  with 
aluminum  oxide  because  of  its  ready  availability.  The  experiments  were  performed  using  an 
unmodified  optical  microprobe^  to  both  excite  the  luminescence  and  to  collect  and  analyze 
the  resulting  fluorescence  spectrum  using  an  attached  spectrometer.  (The  optics  of  the  Mi¬ 
croprobe  have  been  described  in  an  earlier  publication  [9]).  A  number  of  laser  frequencies 
have  been  used  to  excite  the  luminescence  but  the  results  reported  here  were  all  obtained 
using  an  argon  ion  laser  with  a  wavelength  of  514  nm.  The  experimental  procedure  was  to 
select  a  region  of  interest  on  the  sample  with  the  attached  optical  microscope,  focus  the  laser 
beam  to  a  spot  on  the  selected  feature  and  analyze  the  fluorescence  produced.  The  intensity 
of  the  Ri  and  Ra  fluorescence  lines  were  typically  scanned  by  integrating  over  0.5  second 
intervals  at  a  spacing  of  0.2  wavenumbers,  with  the  intensities  being  recorded  under  com¬ 
puter  control.  The  collected  data  were  subsequently  analyzed  with  curve  fitting  algorithms, 
included  in  the  LabCalc  software  package^,  to  identify  the  position  of  the  fluorescence  peak. 
By  using  objective  lenses  of  50  x  and  100  x  magnifying  powers  minimum  spot  sizes  of  ~  5/im 

‘  tiutruiiMntt  SA,  FUman  Microprobe  Model  T64000. 

^Gaiictic  Indiutriea  Corp. 
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and  ~  I  -2/im  can  be  produced.  Larger  probe  sizes  were  formed  by  partial  defocusing.  .Also, 
by  appropriate  choice  of  the  collection  aperture  size,  the  attainable  axial  resolution  could  be 
varied  by  the  collection  optics  from  ~  lO/im  to  ~  200/im  [10]. 

■All  the  measurements  were  made  with  the  sample  at  room  temperature.  .Although  a 
change  in  temperature  is  known  to  cause  a  line  shift  of  ~  0.14  cm~V“C  [11]  (to  smaller 
wavenumbers  with  increasing  temperature),  no  line  shift  was  noted  when  the  spot  size  was 
systematically  decreased  (thereby  increasing  the  power  density  incident  on  the  sample)  indi¬ 
cating  that  no  significant  heating  of  the  probed  volume  occurred.  This  lack  of  sensitivity  to 
spot  size  was  attributed  to  both  the  low  absorption  of  ruby  at  the  excitation  frequencies  and 
the  relatively  low  illumination  intensities  used.  Variations  were,  however,  noted  with  varia¬ 
tions  in  room  temperature  and  so  corrections  were  made  to  the  peak  shifts.  This  involved 
corrections  for  the  change  in  room  temperature  using  the  above  temperature  dependence  of 
line  frequency  and  changes  of  the  spectrometer  dimensions  as  a  result  of  thermal  expansion. 
The  latter  was  corrected  by  simultaneously  monitoring  a  characteristic  Neon  line  at  14431 
cm“*. 

4,1  Polycrystalline  Sample  Under  Stress 

In  order  to  test  the  validity  of  equations  7  and  8  concerning  the  stress  dependence  of  the 
peak  shift  and  broadening  in  a  polycrystalline  material,  polished  bars  of  polycrystalline  alu¬ 
mina  (Coors  AD96)  were  subject  to  an  applied  stress  and  the  fluorescence  spectra  recorded. 
The  stress  was  applied  using  a  four-point  bend  stage  specially  designed  to  fit  under  the  op¬ 
tical  microscope  of  the  microprobe.  The  polycrystalline  alumina  contained  sufficient  levels 
of  chromium  impurity  to  provide  sufficiently  strong  fluorescence  peaks  to  be  analyzed.  In 
order  to  average  over  many  grains,  a  probe  diameter  of  ~  50/im  was  used  for  the  fluorescence 
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measurements.  The  grain  size  of  the  alumina  was  2-20/rm. 

Since  alumina  is  thermally  anisotropic,  the  analysis  in  section  III.  1  has  to  be  extended  to 
separate  the  effect  of  thermal  residual  stress  existing  in  the  sample  from  any  applied  stress. 
The  shift  from  a  small  volume  element  dv  (much  smaller  than  the  grain  size)  is 

Ai/  =  2n,{a^  +  al)  +  n^icT-^  +  aj)  (24) 

where  the  superscripts  ,4  and  T  denote  applied  and  thermal  stresses  respectively.  Since  the 
above  equation  is  linear,  the  effects  of  two  types  of  stress  can  be  separated: 

Ai/  =  Ar/-^  +  Ai/^  (25) 


Therefore,  the  mean  shift; 


At'  = 


+  Ai/^)di;  =  At''^  +  At'^ 


(26) 


Since  the  shift  generated  by  thermal  stress  (see  section  4  3)  can  be  measured  independently. 
At'-’’  can  be  obtained. 

The  broadening  is  calculated  by  using  equation  .A8: 


=  :^y(Ai'^  +  Az/^)^di;- Ai/' 

=  <  (Ai''^)^  >  +  <  (Ai/^)^  >  +2  Aiz-^At/^  -  At'''  Ai/^l 


(27) 


So  the  broadening  can  be  separated  provided  that  At'"^  and  At/^  are  uncorrelated. 

The  measured  peak  shifts  for  the  Ri  and  R?  fluorescence  lines  as  a  function  of  applied 
stress  are  shown  in  figure  2.  The  lines  through  the  data,  fitted  by  a  least  square  analysis. 
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shows  that  the  stress  dependence  of  the  shift  of  the  Ri  and  R2  fluorescence  lines  is  2.46  and 
2.50  cm“'GPa~‘  respectively.  The  flt  coefficient  R  is  0.997  for  both  lines.  For  comparison 
to  the  theory  in  section  III.l,  equation  7  would  predict  that  when  a  uniform  stress  is  applied 
to  polycrystalline  alumina,  the  peak  shift  would  be  given  by: 

A^  =  ^(2n,  +  nj(T^ 

Using  the  hydrostatic  pressure  dependence  of  peeik  shifts  for  the  Ri  and  Rj  lines  given  by 
Munro  et*  [12],  the  sum  of  the  spectroscopic  coefficients,  2na  +  He  have  the  values  of 
7.59  and  7.62  cm~‘  GPa~*  for  the  Ri  and  R2  lines  respectively.  The  peak  shifts  of  the 
Ri  and  Rj  lines  would  then  be  2.53  and  2.54  cm“’GPa~\  close  to  that  observed.  The 
uniaxial  compression  experiment  by  Feher  and  Sturge  [13]  provided  values  for  Ho  and  Oc 
independently  with  values  of  2.7  and  2.15  cm"^  GPa~^  for  both  Ri  and  Rj  lines.  Using 
these  data,  the  predicted  stress  dependence  would  be  2.52  cm“'GPa“‘,  again  close  to  that 
observed.  VV'e  can  therefore  conclude  that  our  measured  shifts  fall  within  the  uncertainty  of 
the  existing  data  on  the  piezo-spectroscopic  coefficients. 

The  measured  peak  widths  for  the  Rx  and  R2  fluorescence  lines  as  functions  of  applied 
stress  are  shown  in  figure  3.  For  this  material  the  systematic  broadening  due  to  the  applied 
stress  of  both  lines  is  small  compared  to  the  variability  from  one  point  on  the  sample  to 
another. 

4.2  Single  Crystal  With  Surface  Stress 

\  single  crystal  treated  to  produce  a  surface  stress  represents  a  special  case  of  the  inhomo- 
geneously  stressed  single  crystal  considered  in  section  III. 2.  To  provide  an  experimental  test 
of  the  effect  of  inhomogeneous  stressing,  a  residual  surface  stress  was  produced  in  a  single 
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crystal  of  c-axis  sappbire  by  in-diffusion  of  chromium  and  the  fluorescence  measured  as  a 
function  of  an  applied  stress.  Since  chromium  oxide,  Cr^Oa,  is  isostructural  with  .A.I2O3  but 
has  a  larger  lattice  parameter,  diffusion  of  chromium  into  the  surface  of  AI2O3  is  expected 
to  produce  a  compressive  surface  stress. 

As  chromium  is  both  the  fluorescing  ion  and  the  ion  causing  the  surface  stress,  the  analysis 
of  section  111.2  must  include  the  stress  created  by  the  in-diffusion  of  chromium.  In  this  case, 
the  frequency  shift  is  expected  to  be  the  sum  of  three  distinct  contributions,  that  due  to 
axiy  applied  stress,  that  due  to  diffusion-induced  stresses  <r^  and  that  due  to  variations  in 
fluorescing  ion  concentrations: 

At/  =  n.2  crtj  +  n.2  trj  +  J(c  -  co)  (28) 

The  overall  peak  shift  from  the  in-diffused  layer  can  be  written  as: 

Ai/ =  y  Ai/ c(2)^(2)<iz  =  At/'"' -f  Ai/^ -h  At/*"  (29) 

Since  the  frequency  shift  in  equation  28  is  the  sum  of  the  individual  components,  the  mean 
peak  shifts  of  each  component  are  also  linearly  independent  as  shown  by  equation  29.  Thus, 
in  an  experiment  in  which  an  external  stress  is  applied  to  a  residually  stressed  surface,  the 
mean  frequency  shift  observed  is  predicted  to  be  due  to  the  applied  stress  only; 

A:/^  =  j  c{z)g{z)dz  =  n.j  cr.'J  (30) 

For  four-point  bending,  the  applied  stress  is  essentially  constant  over  the  depth  of  the 
in-diffused  layer,  and  so  the  mean  frequency  shift  will  be  linearly  related  to  the  applied 
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stress. 


For  the  experiment,  iu-diffusion  of  chromium  was  performed  in  a  three  step  process.  In 
the  first  step,  a  chromium  acetate  solution  was  spun-coated  onto  a  c-axis  oriented  sapphire 
disc.  The  coating  was  then  pyrolyzed  and  oxidized  by  heating  to  600'’C  in  air.  Then, 
finally,  the  disc  was  annealed  at  IfiOO’C  to  diffuse  in  the  chromium  and  beams  cut  from  the 
disc.  To  investigate  the  effect  of  applied  stress  the  sapphire  beam  was  loaded  in  four- point 
bending  with  c-axis  perpendicular  to  the  loading  plane,  placed  under  the  microscope  of  the 
microprobe,  and  the  resulting  fluorescence  recorded  as  a  function  of  load.  The  observed  peak 
shift  for  the  Ri  and  R2  lines  as  a  function  of  applied  stress  is  shown  in  figure  4.  For  these 
data,  the  stress  dependence  is  2.42  and  2.61  cm~‘GPa"‘  for  the  Rj  and  R2  lines  respectively, 
which  compares  favorably  to  the  mean  Ha  value  of  2.7  cm^'GPa”'  given  in  [13].  The  fit 
coefficient  R  equals  0.996  and  0.994,  respectively. 

4.3  Thermal  Expansion  Anisotropy  Stresses 

Residual  stresses  are  known  to  be  produced  in  polycrystalline  alumina  on  cooling  as  a 
result  of  the  crystallographic  anisotropy  in  thermal  expansion  coefficients  along  the  c-  and 
a-  axes.  In  order  to  measure  the  average  residual  stresses,  the  fluorescence  from  three 
different  sintered,  polycrystalline  aluminM  (Coors  .ADOS,  Coors  .^D  995  and  LCSBl)  and  a 
polycrystalline  sJumina  hipped  at  UCSB  (UCSB2)  was  recorded  using  an  excitation  probe 
diameter  of  50^m.  The  UCSBl  material  is  a  high  purity  alumina  sintered  at  ISOO^C  for  1  hr. 
in  hydrogen.  The  probe  size  was  much  larger  than  the  grain  size  of  the  ceramics,  assuring 
that  the  fluorescence  was  being  averaged  over  a  large  number  of  grains.  1  he  ceramic  samples 
contained  an  adequate  chromium  impurity  level  for  the  Ri  and  R2  fluorescence  peaks  to  be 
recorded  with  sufficient  signal-to-background  ratio  that  precise  measurements  could  be  made 
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of  the  peak  shape  and  peak  shift.  The  peak  shift  and  peak  broadening  were  measured  relative 
to  those  recorded  from  a  polished  single  crystal  of  sapphire  that  also  contained  a  trace  level 
of  chromium  impurity.  For  comparison,  the  spectrum  obtained  from  the  .^0995  sample  is 
plotted  together  with  that  of  a  single  crystal  sapphire  in  Figure  5.  Several  spectra  were  taken 
at  different  spots  for  each  sample,  the  average  peak  shifts  and  broadening  for  the  samples 
are  listed  in  table  I.  Also  shown  in  the  table  are  the  average  stresses,  aZ,  Wl  and  the  breath 
of  stress  distribution  (Tu,  calculated  using  equations  16,  17  and  20  above. 

5  CLOSING  REMARKS 

The  preceding  analyses  and  experiments  illustrate  how  the  technique  of  fluorescence  spec¬ 
troscopy  can  be  used  to  make  measurements  of  commonly  arising  types  of  .stresses  in  cerzim- 
ics.  This  is  possible,  despite  the  fact  that  the  fluorescence  shift  is  a  scalar  quantity  whereas 
the  complete  stress  state  is  tensorial,  because  considerable  simplification  of  the  equations 
(equation  .3)  can  be  made  in  many  practical  situations.  For  instance,  stresses  can  be  applied 
at  prescribed  orientations  to  the  crystallographic  axes  of  single  crystals.  Also,  the  poly¬ 
crystalline  nature  of  many  ceramics  enables  the  effects  of  crystallographic  orientation  to  be 
averaged  out.  For  instance,  as  equation  7  demonstrates,  the  frequency  shift  from  a  stressed 
polycrystalline  ceramic  is  independent  of  its  crystal  structure  and  is  only  dependent  on  the 
trace  of  me  applied  stress. 

Two  of  the  analytical  results  are  of  particular  interest.  One,  is  the  finding  that  the 
fluorescence  shift  is  dependent  only  on  the  applied  stress  and  not  on  the  residual  stresses 
produced  by  either  thermal  expansion  mismatch  as  in  the  case  of  polycrystalline  AI2O3,  or 
by  ion  in-diffusion  as  in  the  case  of  Cr^'*'  in-diffusion  into  sapphire  single  crystal  .  The  later 
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also  suggests  that  a  fluorescent  dopant  ion  can  be  deliberately  introduced  into  a  structure  or 
material  to  monitor  the  applied  stress  through  changes  in  the  fluorescence  frequency.  In  fact, 
we  have  been  utilizing  this  characteristics  to  make  measurements  of  stresses  in  metal-ceramic 
laminates  [14]. 

The  second  notable  result  concerns  the  measurement  of  the  residual  stresses  due  to 
anisotropic  thermal  expansion.  The  analysis  indicates  that  the  mean  stress  in  the  a-  and 
c-  crystallographic  directions  in  aluminum  oxide  can  be  obtained  from  the  mean  frequency 
shift,  and  that  the  distribution  in  internal  stresses  can  also  be  determined.  These  results  rep¬ 
resent  an  advance  on  previous  studies  of  anisotropic  thermal  expansion  stresses  in  aluminum 
oxide  [2, 15).  As  already  mentioned,  Grabner's  work  led  to  an  incomplete  expression  for  the 
relationship  between  the  frequency  shift  and  stress,  and  so  the  values  he  derived  are  probably 
unreliable.  The  work  of  Blendell  and  Coble  [I5j,  was  primarily  concerned  with  the  diffusional 
relaxation  of  the  anisotropic  thermal  expansion  stresses.  They  assumed  an  average  stress  in 
the  material  and  calculated  its  relaxation  for  different  grain  sizes  and  cooling  rates.  They 
did  not  relate  this  average  stress  to  any  of  the  crystallographic  stresses  nor  did  they  describe 
how  they  determined  this  average  stress  from  their  fluorescence  observations.  However,  the 
measurements  in  table  I  indicate  that  the  mean  stresses  in  the  samples  of  polycrystalline 
alumina  examined  axe  similar  in  magnitude  to  those  calculated. 
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APPENDIX 


A  Peak  Shifts  and  Broadening  due  to  Stresses 

If  the  probed  volume  were  strain  free,  the  fluorescence  peaks  would  be  Lorent^ian  in  shape 
and  have  a  breadth  dependent  only  on  the  finite  temperature  of  the  sample.  The  presence 
of  stresses,  however,  causes  the  peaks  to  be  shifted  and  dispersed  as  shown  in  the  following. 

Let  the  peak  shape  in  the  unstrained  crystal  be  represented  by  the  frequency  function, 
f{v),  with  the  peak  centroid  given  by: 


f{v)udu 


(Al) 


and  the  area  of  the  peak  is  normalized  to  be  unity: 


(A2) 


The  dispersion  of  the  peak,  corresponding  to  the  peak  width,  can  be  expressed  by  the 
second  moment  of  the  distribution: 


<  >0=  j (A.3) 

The  effect  of  any  stress  distribution  can  be  related  to  measured  peak  shape  by  assuming 
that  the  probability  of  causing  a  frequency  shift,  Ai/,  can  be  expressed  as  a  weighting  function 
The  resulting  peak  is  then  given  by  the  function: 


F(u)  =  J  /(i/- 


(A4) 
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Since  /  =  /  F{i/)du  =  1  the  area  under  the  peak  is  conserved. 

The  centroid,  or  mean  value,  of  the  peak  from  the  stressed  volume  is  given  by  the  first 
moment  of  the  function  representing  the  peak; 

Uc  =  j  F{u)t/du  {.\o) 

The  peak  shift  is  then: 

At/  =  —  J/Q 

=  j  (.^6) 

As  spectra  f{u)  and  F{u)  can  be  obtained  from  unstressed  and  stressed  regions,  the  mean 

shift  ^1/  can  be  determined  and  from  that  the  stress  deconvoluted. 

The  width  of  the  stress  broadened  peak,  F{v),  is: 

<  1/2  >p=  j  F{i'){v  -  uc)^du  (A7) 

The  broadening  is  then  seen  to  be  the  difference  between  the  measured  widths  of  the  peaks 
from  the  strain  free  regions  and  from  the  stressed  regions: 

<  >F  —  <  >0  =  j  W(J:^Li/)dAu  ~ 

=  <  Ai/^  >  (A8) 

Again,  the  measured  broadening  and  the  stress  distribution  is  related  by  the  weighting 
function. 
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B  Polycrystalline  Averages  of  Direction  Cosines 

Collected  here  for  convenience  are  the  polycrystalline  averages  of  the  direction  cosines  a 
appearing  in  tensor  transformation  relations  of  the  form; 


'Bi: 


In  terms  of  the  usual  Euler  angles  (d,  (^,  tp),  the  transformation  matrix  is 


a„  = 


cos  <f>  cos  0  —  sin  <p  cos  d  sin  0 
sin  4>  cos  0  +  cos  0  cos  d  sin  0 
sin  ^  sin  0 


cos  0  sin  0  —  sin  0  cos  6  cos  0 
sin  0  sin  0  +  cos  0  cos  $  cos  0 
sin  B  cos  0 


sin  d  sin  0 
~sinflcos0 
cos  9 


where 


O<0<27r,  0<^<?r,  O<0<2:r 


(B2) 


If,  as  in  a  polycrystalline  material,  the  orientations  of  the  grains  are  at  random,  the 
probability  a  grain  director  is  oriented  at  (^,  0,  0)  is  given  by 


0, 0)  =  ■^-rs\n9d9d<t)dTp 

OK* 


(B3) 


By  averaging  over  all  possible  angles,  the  following  relations  used  to  derive  equation  7 
and  8  are  obtained: 


nqflifc  —  0  ,  j  ^  k 
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<^21^22^23  —  (22ia22<223  —  021022^23  —  0 


22^^31^32  =  <l2l023<^3l033  =  <^22023*^32^33  = 
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figure  captions 

1.  Schematic  diagram  of  a  representative  grain  cut  by  an  arbitrary  dividing  plane  (c  =  0) 
through  a  polycrystalline  material,  illustrating  the  coordinates  and  angles  used  in  the  text 
in  describing  the  orientation  of  the  grain. 

2.  Measured  peak  shifts  of  the  Ri  (top)  and  R2  (bottom)  Buorescence  lines  from  a 
polycrystalline  alumina  sample  as  a  function  of  applied  stress.  The  sample  was  loaded 
under  four-point  bending  and  spectra  were  recorded  from  different  positions  on  both  the 
tensile  and  compressional  sides  and  over  the  stress  range  indicated.  The  solid  line  through 
the  data  points  represents  the  line  of  best-fit.  Its  slope  is,  within  experimental  uncertainty, 
the  same  as  that  predicted  using  equation  7. 

.3.  Measured  peak  broadening  of  the  Ri  and  R2  fluorescence  lines  as  a  function  of  stress. 
The  broadening  was  obtained  from  the  same  spectra  as  were  used  to  make  the 
measurements  in  figure  2. 

4.  Peak  shift  of  the  Rj  and  R2  fluorescence  lines  as  a  function  of  applied  stress  for  a 
sapphire  crystal  containing  a  surface  stressed  layer.  The  solid  line  is  the  line  of  best  fit  to 
the  data. The  surface  stress  was  produced  by  in-diffusion  of  chromium.  Details  are  given  in 
the  text. 

5.  .\  comparison  of  the  Ri  and  R2  fluorescence  lines  from  sapphire  and  a  polycrystalline 
alumina  (Coors  AD995).  The  broadening  of  the  lines  from  the  polycrystalline  alumina  is  a 
result  of  the  distribution  of  residual  thermal  expansion  anisotropy  stresses  within  the 
material.  Both  materials  contain  a  trace  of  chromium  impurity,  sufficient  to  cause 
fluorescence. 
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TABLE  I 


SAMPLE 

3:17 

cm"' 

<ya 

MPa 

57 

.MPa 

MPa 

>/<  Aj/2  > 

cm"' 

AD96 

-0.21 

-190 

380 

2.5 

430 

AD995 

-0.15 

-135 

270 

3.6 

620 

UCSBl 

-0.15 

-135 

270 

1.9 

330 

UCSB2 

-0.23 

-210 

420 

1.4 

240 
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Abstract 

The  residual  stresses  in  c*axis  sapphire  fibers  in  a  Y  -TiAl  matrix  and  in  a 
polycrystaUine  AI2O3  matrix  as  a  function  of  distance  below  a  surface  are 
determined.  They  are  obtained  from  the  shift  in  frequency  of  the  characteristic 
R2  fluorecence  line  of  chromiiun  in  sapphire  obtained  by  foccusing  an  optical 
probe  at  different  depths  in  a  sapphire  fiber  intersecting  the  surface  of  the 
composite.  The  metkiod  is  described  together  with  its  calibration.  Both  the  axial 
and  radial  components  of  the  residual  stress  in  the  fiber  are  observed  to  vary 
over  a  depth  of  approximately  the  fiber  diameter  and  are  then  almost 
independent  of  depth. 


1  Introduction 


When  a  composite  material  is  cooled  down  from  its  fabrication  temperature,  therrr^al 
mismatch  residual  stresses  develop,  which  can  have  a  profound  influence  on  the  mechanical 
properties  of  the  composite.  For  example,  in  many  ceramic  fiber  reinforced  metal  matn.x  or 
glass  matri.x  composites  the  thermal  expansion  coefficient  of  the  matrix  is  considerably 
larger  than  that  of  the  fiber  O/  and  so  large  compressive  stresses  can  develop  in  the  fiber. 
The  axial  compressive  stress  creates  a  large  stress  difference  between  debonded  and  bonded 
regions  of  the  interface  when  the  material  is  under  applied  tension,  facilitating  interfacial 
crack  growth  [l,2|.  a  pre-requisite  for  fiber  pull-out.  On  the  other  hand,  the  large  radial 
compressive  stress  produces  large  friction  forces  in  'he  debonded  region,  which  retards  fiber 
sliding  and  consequently  shortens  fiber  pull-out  length  [2,3].  U  is  therefore  possible  to  opti¬ 
mize  the  mechanical  properties  of  such  composites  by  tailoring  the  residual  stresses  through 
appropriate  processing  conditions  and  the  use  of  interfacial  coatings  [4]. 

Unfortunately,  accurate  information  about  the  residual  stresses  in  embedded  fibers  is  usu¬ 
ally  unavailable  since  their  measurement  is  particularly  difficult.  .\s  a  result,  estimations 
based  on  the  thermal  and  elastic  properties  of  the  individual  material  components  are  usu¬ 
ally  used  in  various  analyses.  This  approach  requires  accurate  knowledge  of  the  effective 
temperature  change  AT,  which  is  usually  an  unmeasurable  quantity,  and  is  difficult  to  es¬ 
timate  without  detailed  knowledge  of  the  stress  relaxation  mechanisms,  especially  near  the 
interface.  Furthermore,  such  estimates  can  not  take  into  account  the  variations  from  fiber 
to  fiber  due  to  irregularities  in  coatings  and  interfacial  roughness:  details  which  are  impor¬ 
tant  in  understanding  individual  fiber  puli-out/push-out  tests.  Direct  measurement  of  the 
residual  stresses  in  fibers  is.  in  principle,  possible  by  a  variety  of  diffraction  techniques.  For 
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instance.  X-  ray  diffraction  using  uitra-fine  beams  can  provide  strain  information  of  individ¬ 
ual  fibers.  However,  it  suffers  from  having  a  relatively  small  penetration  depth  (~  )0  ^.m) 
[5.61  and  can  therefore  only  measure  stresses  near  the  surface,  which  could  be  significantly 
different  from  those  deeper  inside  the  composite.  .Veutron  scattering,  whilst  capable  of  prob¬ 
ing  fiber  stresses  throughout  the  composite  :7,Sj,  has  insufficient  spatial  resolution  to  study 
individual  fibers  or  to  ascertain  stress  distributions  along  the  fibers. 

.An  alternative  method  of  measuring  stresses  in  fibers  is  by  using  piezo-spectroscopic 
effects.  The  basis  of  such  methods  is  that  the  characteristic  luminescence  lines  (Raman 
or  fluorescence)  of  the  material  shift  with  stress.  By  utilizing  an  optical  microscope  to 
define  a  probe,  it  is  relatively  easy  to  select  regions  of  interest  with  spatial  resolution  of 
only  a  few  microns  [9]  either  on  the  surface  of  a  composite,  or  if  the  matrix  is  transparent, 
within  the  composite.  For  example,  in  the  case  of  polydiacetylene  fibers  embedded  in  an 
optically  transparent  epoxy  [lOj.  Raman  spectra  were  collected  from  different  positions  along 
a  fiber  to  obtain  the  stress  distribution  along  it.  Since  optical  spectrum  can  be  collected  and 
analyzed  with  relative  ease  compared  to  most  diffraction  methods,  such  measurements  can 
be  performed  routinely. 

In  this  paper,  we  present  a  new  method  of  measuring  residual  stresses  in  embedded 
transparent  fibers,  one  that  utilizes  optical  fluorescence,  and  illustrate  its  application  to 
determining  the  residual  stresses  in  sapphire  fibers  in  7-Ti.Al  and  .AI2O3  matrices.  The 
theoretical  basis  of  the  fluorescence  method  was  described  in  detail  in  an  earlier  work  [11]. 
Briefly,  the  frequency  shift  of  a  characteristic  fluorescence  line  is  stress  dependent  and  is 
given  by: 

Ai/  =  n,;  cr’  ( 1 ) 
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where  are  the  piezo-spectroscopic  coefficients  relating  frequency  shift  to  stress,  and  7', 
is  the  stress  tensor  represented  in  the  crystallographic  frame  of  reference. 

The  method  is  to  use  an  optical  microscope  to  focus  a  light  beam  on  the  end  of  an 
e.xposed  sapphire  fiber  (figure  1 ).  measure  the  frequency  shift  of  a  characteristic  luminescence 
line  and  then  repeat  the  measurement  at  successively  deeper  points  of  focus  into  the  fiber. 
Using  the  lenses  of  the  microscope  to  form  an  optical  probe  enables  the  luminescence  to 
be  collected,  for  any  particular  point  of  focus,  from  a  relatively  small  depth  of  field.  This 
facilitates  high  spatial  resolution  [L2.l3j  along  the  fiber  but  at  the  same  time  means  that  the 
observed  frequency  shifts  have  to  be  deconvoluted  using  a  previously  established  calibration 
of  the  microscope  convolution  function  (the  depth  of  field  function).  This  is  described  in  the 
following  section.  .Also  described  is  a  test  of  the  accuracy  of  the  stress  profiling  method  in 
w’hich  the  stress  through  the  thickness  of  a  ruby  beam  under  four-point  loading  is  measured 
and  compared  with  the  known  elastic  solution. 

To  illustrate  the  application  of  the  technique  we  have  meaisured  the  residual  stresses 
in  Saphikon  sapphire  fibers.  They  have  sufficient  levels  of  chromium  impurity  to  produce 
sufficiently  intense  fluorescence  lines  Ri  and  Rj  of  ruby  when  e.'ccited  by  an  argon  ion  laser. 
Because  the  fibers  are  single  crystal  with  their  c-axis  along  the  fiber  a.xis.  the  residual  stresses 
inside  the  fibers  are  expected  to  have  cylindrical  sym,metry.  .As  a  result,  there  are  only  two 
independent  stress  components,  the  radial  component  7r  and  the  axial  component  7j.  This 
simplifies  equation  1  to  be: 

Ai/  =  2113  -r  rlr  7;  (  21 

where  fl,  and  flc  are  the  piezo-spectroscopy  coefficients  for  the  a  and  c  directions.  They 
have  been  measured  to  have  values  of  2.70  and  2.15  cm~'GPa~’ ,  respectively  [14].  .A  positive 
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shift  implies  tension  and  a  negative  one  implies  compression.  In  section  .3.  measurements  of 
line  shifts  along  fibers  in  both  7-TiAl  and  .-MaOs  matrices  as  a  function  of  depth  are  made. 
Since  the  frequency  shifts  are  dependent  on  both  the  a.xial  and  radial  stress  components, 
their  separation  requires  the  solution  of  the  appropriate  elastic  equations.  This  is  described, 
using  the  results  of  a  finite  element  computation  in  section  4,  From  this,  the  stresses  as  a 
function  of  depth  are  determined. 

2  The  Through-Focus  Depth  Profiling  Method  ^ 

2.1  Depth  of  Field  Function 

The  depth  resolution  of  the  optical  microprobe  used  for  the  fluorescence  measurements 
is  determined  by  the  numerical  aperture  of  the  objective  lens  and  the  size  of  the  collection 
aperture  used  [12,13].  As  illustrated  in  figure  2.  when  the  excitation  laser  is  focused  at  a 
distance  -  below  the  top  surface  of  a  transparent  material  with  uniform  concentration  of 
fluorescing  species,  only  the  signal  generated  in  a  small  volume  above  and  below  the  focal 
plane  is  collected.  The  total  signal  collected  is  then  the  integral  of  the  fluorescence  from 
different  depths: 

I{z)  s=  j  ^{u)du  (3) 

where  g{u)  is  the  depth  of  field  function  of  the  lens  used  to  describe  the  relative  collection 
efficiency  as  a  function  of  u,  the  distance  from  the  focal  plane. 

By  moving  the  focal  plane  from  a  position  above  the  top  surface  to  a  position  below 
the  surface  of  a  thick  sapphire  disk  {t  —*  oc)  containing  adequate  chromium  concentration. 

'  The  generaJ  procedures  for  opcicai  duorescence  measurement  and  analysis  were  described  in  detail  in  [1  Ij  and  therefore  not 
included  in  this  work. 
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and  recording  the  fluorescence  line  intensity  as  a  function  of  depth  c,  the  depth  resolution 
function  g  is  obtained  from  equation  3,  namely; 


i-i) 


The  depth  resolution  function  for  our  optical  microprobe  using  a  50/0.55  objective  lens 
and  a  200  /im  collection  aperture  was  measured  using  the  above  method  and  is  shown  in 
figure  3.  As  can  be  seen,  the  collection  efficiency  is  maximum  at  the  focal  plane,  and  decrease 
rapidly  with  distance  on  either  side  of  the  focal  plane. 

2.2  Measurement  of  the  Frequency  Shift  Depth  Profile 

VVe  now' suppose  that  there  exists  a  stress  field  that  only  varies  with  depth  z.  From 
equation  2.  the  shift  of  a  fluorescence  line  associated  with  the  increment  in  signal  from  a 
narrow  slice  of  material  at  depth  r  +  u  is 


+•  u)  =  2U.a<^r{^  -h  u)  -r  ric  cr.(r  u) 


io) 


The  measured  shift  with  focal  point  positioned  at  depth  r  is  then  a  weighted  average  of  the 


signals  from  ail  depths; 


— —  fl/  u)giu)du 

^u{z)  — 

9{n)du 


Provided  the  depth  of  field  function  g{z)  is  known,  the  real  depth  profile  Ai/(z)  can  be 
deconvoluted  from  equation  6. 

To  establish  that  the  foregoing  methodology  is  correct,  the  stress  field  through  the  thick¬ 
ness  of  a  dilute  ruby  beam  in  4-poini  loading  was  measured  and  compared  with  the  elastic 
solution  for  the  stress.  .\  number  of  fluorescence  spectra  were  taken  as  the  laser  beam  was 


focused  at  different  depths  from  the  top,  tension  surface  to  the  bottom,  compression  surface 
of  the  beam.  The  measured  shifts  are  plotted  in  figure  4  as  the  solid  dots.  The  deconvoluted 
data  is  shown  by  the  solid  line.  The  deconvolution  was  obtained  by  assuming  that  the  real 
shift  Ai/(r)  is  a  linear  function  of  depth  a  bz  since  it  is  known  that  the  stress  distribution 
varies  linearly  with  the  depth  r  for  a  rectangular  beam  under  4-point  loading.  This  function 
is  then  convoluted  with  the  depth  of  field  function  g(z).  By  fitting  the  resulting  convoluted 
function  to  the  measured  shifts,  the  coefficients  a  and  b  are  obtained.  The  initially  assumed 
function,  with  the  coefficients  obtained  by  fitting,  is  then  regarded  as  the  deconvoluted  shift 
At'(c)  directly  related  to  the  stress  distributions  by  equation  2.  The  straight  solid  line  in 
figure  4  therefore  represents  the  real  shifts.  As  e.xpected,  and  required,  the  zero  shift  po¬ 
sition  coincides  with  the  neutral  plane  of  the  beam.  The  stress  distribution  calculated  by 
using  equation  2  is  also  in  good  agreement  with  that  obtained  from  beam  theory  and  loading 
conditions.  As  a  further  check,  the  convoluted  shift,  obtained  by  convoluting  the  derived 
values  of  the  real  shift  with  the  instrument  function  g{z)  is  shown  by  the  dashed  curve.  It 
has  e.xcellent  fit  with  the  measured  data  except  for  the  slight  deviation  near  the  compression 
surface  which,  we  believe,  is  due  to  the  fact  that  the  focus  deteriorates  after  light  has  passed 
all  the  way  through  the  ruby  beam.  These  observations  thus  confirm  that  the  through-focus 
method  can  indeed  accurately  measure  the  depth  profile  of  the  stresses. 

3  Measured  Frequency  Shifts  From  Embedded  Sapphire  Fibers 

The  depth  profiling  method  has  been  applied  to  two  different  sapphire  fiber  composites: 
the  first,  a  y-Ti.Al  matrix  composite  with  carbon  black/alumina  double  fiber  coatings  on 
the  sapphire  fibers  and,  the  second,  a  polycrystalline  .\l2O3  matrix  with  a  thin  molybdenum 
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coating  on  the  fibers.  In  each  case,  the  fiber  diameter  is  about  120  Thin  slices  of 
composite  with  parallel  surfaces  perpendicular  to  the  sapphire  fibers  were  cut  to  about  loO 
and  600  thick  for  the  *'-TiAl  matrix  and  the  .AIjOj  matrix  composites  respectiveiy. 
The  top  surface  of  the  specimens  were  then  polished  to  obtain  an  optical  finish  at  the  fiber 
ends. 

The  measured  and  the  deconvoluted  frequency  shift  as  a  function  of  depth  into  a  sapphire 
fiber  in  7-Ti.-\l  matrix,  as  well  as  the  convoluted  function  are  plotted  in  figure  5.  Since 
7-Ti.\l  has  considerably  larger  thermal  expansion  coefficient  than  sapphire  {table  1).  the 
residual  stress  developed  by  thermal  mismatch  in  the  fiber  should  be  compressive  in  both 
radial  and  axial  directions.  This  is  consistent  with  the  large  negative  shifts  observed  when 
the  probe  is  focused  deep  in  the  interior  of  the  specimen.  .Vear  the  surface,  the  axial  stress 
approaches  zero  as  required  by  force  balance.  The  deconvoluted  shift  becomes  positive 
near  the  surface  indicating  that  the  radial  stress  becomes  tensile  near  the  surface.  In  the 
case  of  the  sapphire  fibers  in  polycrystalline  .\l2O3,  the  residual  stress  is  mainly  caused 
by  the  anisotropy  of  the  thermal  expansion  coefficient  of  sapphire.  .As  show  in  figure  6,  the 
relatively  large  negative  shifts  near  the  surface  indicate  that  the  radial  stress  near  the  surface 
is  compressive,  since  the  axial  stress  must  approach  zero  due  to  force  balance  requirement, 
and  thus  does  not  contribute  to  the  shift.  The  magnitude  of  the  negative  shift  becomes 
progressively  smaller  with  increasing  of  depth  implying  that  the  negative  shift  associated 
with  the  radial  compressive  stress  is  partially  compensated  by  the  positive  shift  resulting 
from  the  build  up  of  the  axial  tensile  stress.  These  observations  are  consistent  with  the  fact 
that  the  thermal  expansion  coefficient  of  polycrysialline  .\l2O3  is  larger  than  that  of  the 
sapphire  in  the  a-direction.  but  smaller  than  that  in  the  c-direction.  .A  notable  feature  of 
the  depth  profile  of  the  shift,  in  both  systems,  is  the  near  exponential  dependence.  This 
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is  inciicati%'e  of  the  fibers  being  fully  bonded  to  the  matrices  and  the  fiber/matrix  systems 
being  elastic. 

4  Analysis 

The  residual  stresses  in  the  sapphire  fibers  are  related  to  the  observed  frequency  shifts 
presented  in  the  previous  section  by  equation  2.  In  order  to  obtain  actual  values  of  the 
two  stress  components,  an  additional  relationship  between  the  two  components  is  evidently 
required.  Two  difficulties  are  associated  with  establishing  such  a  relationship.  First,  the 
residual  stresses  originate  from  thermal  expansion  mismatch  but  the  temperature  on  cooling 
at  which  stress  relaxation  no  longer  is  effective,  and  stresses  begin  to  build  up,  is  not  known. 
Secondly,  existing  equations  relating  the  residual  stresses  in  a  fiber  embedded  in  an  infinite 
matrix  [2]  to  the  thermal  mismatch  strain  :^a^T  are  unfortunately  not  valid  for  the  thin 
specimens  used  in  this  work.  Furthermore,  rigorous  analytical  solutions  for  the  residual  stress 
field  are  not  available  and  the  approximate  models  used  in  analyzing  stresses  in  fibers,  e.g. 
those  in  analyzing  fiber  pull-out/push-out  tests  ^lo-lT]  breakdown  near  the  fiber  end.  We 
have  therefore  had  to  resort  to  finite  element  calculations  of  the  residual  stress  field  for  the 
sapphire/7-TiAl  and  sapphire/ AlaOa  systems  and  use  an  iterative,  self  consistent  procedure 
to  calculate  the  residuad  stresses  and  the  temperature  at  which  the  stresses  become  frozen 
in.  This  is  described  in  the  following  paragraphs.  First,  however,  it  should  be  noted  that 
the  stresses  can  be  estimated  directly  from  equation  2.  For  the  7-  Ti.Al  matrix  system, 
the  residual  stresses  in  the  fiber  are  expected  to  be  compressive.  If  it  is  assumed  that  the 
magnitude  of  the  radial  component  is  similar  to  the  axial  component  in  the  interior  of  the 


3 


composite,  then  using  the  shift  from  deep  in  the  fiber  we  have  from  equation  2 


.tnt 


2lla  +  He 


~  500  MPa 


For  the  AI2O3  matrix  system,  due  to  the  anisotropy  of  the  thermal  expansion  coefficient  of 
sapphire,  the  fiber  should  be  under  axial  tension  and  radial  compression.  If  we  assume  that 
the  radial  compression  stress  is  about  the  same  in  the  interior  and  near  the  surface,  we  again 
can  estimate  that 


2n, 


~  -100  MPa 


and 


,int 


He 


-  200  MPa 


Quantitative  stress  distributions  can  be  obtained  by  comparing  in  detail  the  measured 
shifts  with  those  calculated  from  the  elastic  field  in  the  fiber/matrix  system.  Since  the  fiber 
coatings  are  relatively  thin  in  both  composites,  we  expect  their  effects  on  the  elastic  field  to 
be  negligible  compared  to  those  of  the  fiber  and  matrix.  Therefore,  in  the  calculations,  the 
elastic  system  is  simplified  as  consisting  solely  of  an  isotropic  matrix  of  modulus  and 
Poisson's  ratio  i/m  and  a  fiber  of  modulus  E/  and  Poisson's  ratio  t/j  (Table  1).  The  specimen 
was  represented  by  a  thin  disk  of  the  same  thickness  as  the  corresponding  specimen  with 
one  cylindrical  fiber  sitting  at  the  symmetry  axis.  The  diameter  of  the  disk  was  chosen  to 
be  about  ten  times  that  of  the  fiber  in  order  to  be  consistent  with  the  low  fiber  volume 
fraction  of  both  composite  systems.  The  thermal  expansion  coefficients  of  the  sapphire  fiber 
in  both  the  radial  direction  0.3  and  axial  direction  Oc  are  also  listed  in  Table  1.  The  mean 


thermal  expansion  coefficient  for  7-TiAl  is  known  to  be  1 1  —  13  x  10~®/“C  depending  on 


the  stress  frozen-in  temperature.  However,  that  of  the  polycrystailine  .\l2O3  is  e.xpected 
to  vary  from  sample  to  sample  depending  on  the  residual  matri.x  porosity.  .\s  a  result  of 
the  complicated  stress  rela.xation  processes  that  occur  during  cooling,  particularly  near  the 
interface,  the  etfective  thermal  expansion  coefficient  ct^  for  both  matrices  and  the  effective 
temperature  change  AT  were  adjusted  in  the  computation  to  obtain  the  best  agreement 
with  the  experimental  frequency  shift  data.  The  following  procedure  was  adopted.  First, 
reasonable  estimates  for  the  matrix  thermal  expansion  coefficient  and  temperature  change 
and  were  used  in  the  finite  element  calculation.  The  resulting  stresses  and 

were  then  used  in  equation  2  to  calculate  the  shift  as  a  function  of  z,  Ai/"j(r).  This 
function  and  the  measured  shifts  are  then  normalized  by  their  values  in  the  interior  of  the 
material:  • 

(7) 

By  comparing  the  calculated  and  measured  normalized  shifts  for  a  series  of  values  of  a“‘, 
the  value  that  gives  the  best  agreement  is  considered  as  the  true  Qrn.  Then,  by  requiring  that 
the  calculated  shift  be  equal  to  the  measured  shift  in  the  interior,  the  effective  temperature 
difference  is  proportion^Llly  scaled: 


AT  =  Ar"‘ 

Ai/i: 


Finally,  the  stress  distributions  are  obtained: 


AT 

Ar>‘  ' 


<^1 


AT 

AP^‘  ■ 


(9) 
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Since  the  thermal  expansion  coefficient  of  7-TiAl  is  considerably  larger  than  a,  and  a,  of 
sapphire,  the  difference  between  and  was  unimportant.  As  a  result,  it  was  observed 
that  the  fitting  of  the  normalized  shift  was  not  sensitive  to  the  actual  value  of  Figure 
7  illustrates  the  fitting  between  the  measured  and  the  calculated  normal  shifts  for  =  12. 
The  effective  temperature  drop  AT  calculated  from  equation  S  is  ~  7.50®C.  which  agrees 
with  a^-T  curve  for  7-Ti.Al  [iSj.  This  temperature  change  is  about  .300®C  smaller  than 
the  difference  between  the  processing  temperature  and  the  room  temperature  indicating 
significant  stress  relaxation  occurring  during  cooling.  The  stress  distributions  along  the 
fiber  calculated  from  equation  9  are  plotted  in  figure  8.  The  axial  stress  approaches  zero  as 
expected  and  the  radial  component  is  positive  at  the  surface  as  was  indicated  by  the  frequency 
shift  measurement.  The  differences  between  the  calculated  and  the  measured  normalized 
shifts  as  shown  in  figure  7  are  due  to  the  simplifications  used  in  the  model,  e.g.,  neglecting 
the  effects  of  coatings  and  eissuming  perfect  interfaciai  bonding.  However,  the  finite  element 
calculations  shows  that  the  stress  distributions  in  the  interior  of  the  specimen  are  not  sensitive 
to  small  changes  in  the  value  of  the  elastic  modulus  and  the  thermal  expansion  coefficient 
of  the  matrix.  This  is  because  the  values  of  the  stresses  are  bounded  by  the  measured 
shifts  through  equations  2  and  scaled  by  equations  8  and  9.  These  observations  suggest 
that  the  stress  distributions  deduced  from  the  fluorescence  data  and  illustrated  in  figure  8 
are  relatively  zwtcurate  especially  in  the  interior  of  the  specimen.  The  possibility  of  some 
interfaciai  debonding  will  be  discussed  in  more  detail  in  the  next  section. 

.\s  expected,  the  normalized  shift  calculated  for  the  .AI2O3  matrix  system  depends  much 
more  sensitively  on  the  value  of  q..;  used  in  the  computation.  .As  show  in  figure  9.  Qm  =  3.7.3 
gives  the  best  fit  between  the  calculated  and  the  measured  normalized  shifts.  The  effective 
temperature  difference  obtained  using  equation  8  is  about  1300®C.  which  is  about  200“C 
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lower  than  the  processing  temperature  again  as  indication  of  stress  relaxation  occurring 
during  cooling.  The  stress  distributions  in  the  sapphire  fiber  are  plotted  in  figure  10.  The 
observation  that  the  calculated  stress  distributions  in  the  interior  are  not  sensitive  to  changes 
of  also  applies  in  this  system.  The  values  of  stress  in  the  interior  differ  less  than  59c  for 
the  three  Om's  shown  in  figure  9.  Since  the  shape  of  the  shift  profile  depends  sensitively  to 
the  change  of  On,  the  fitting  in  figure  9  also  provides  means  of  estimating  the  local  thermal 
expansion  coefficient  of  the  .-\l2O3  matrix,  which  could  be  different  from  the  average  value 
of  the  whole  sample  because  the  density  of  the  matrix  is  slightly  lower  near  the  fibers.  The 
resulting  value  of  can  then  be  used  in  other  analyses,  such  as  for  push-out  tests. 


5  Discussion 


The  analysis  in  the  last  section  is  based  on- the  assumption  that  the  fiber/matrix  systems 
are  elastic  without  any  debonding  occurring  at  the  interface.  This  requires  that  the  inter- 
facial  shear  strength  be  larger  than  the  maximum  interface  shear  stress.  Tc  examine  this 
assumption,  the  interfacial  shear  stress  is  also  calculated  in  the  finite  element  computation. 
The  results  corresponding  to  the  best  fit  conditions  determined  in  the  last  section  are  plotted 
for  the  sapphire/7-TiAl  and  the  sapphire/ AI2O3  systems  in  figure  11.  In  both  systems,  the 
magnitude  of  the  shear  stress  hats  a  maximum  at  the  surface  and  decreases  to  zero  with 
increasing  depth.  This  is  expected  from  the  stress  distributions  in  figures  3  and  10  and  the 
force  balance  requirement: 


r  = 


Rf  da, 

Y~dz 


(10) 


Whilst  the  calculated  interfacial  shear  stress  for  the  sapphire/ .\l2O3  system  is  smaller 
than  the  interfacial  shear  strength  for  the  molybdenum/sapphire  interface  100  —  300 
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MPa)  *l9j  and  therefore  insufficient  to  produce  any  debonding,  th**  ihe’^  .ess  calculated 
for  the  sapphire/7-Ti.AI  system  is  considerably  larger  than  any  reasonable  value  for  the  .ar- 
bon/sapphire  interface  (~  100  —  150  MPa)  [-l.lSj.  It  is  likely  that  ’he  thin  carbon  coating 
on  this  fiber  is  not  continuous  so  that  the  sapphire  fiber  has  direct  contact  with  the  poly- 
crystalline  -AlaOi  outer  coating  to  form  a  strongly  bonded  interface.  This  has  been  oTen 
observed  in  the  same  composite  specimen  during  push-out  tests  [IS].  As  mentioned  in  the 
last  section,  the  fitting  between  the  calculated  and  measured  frequency  shifts  near  the  sur¬ 
face  is  not  perfect.  .As  shown  in  figure  7,  the  calculated  shift  changes  more  rapidly  than 
that  measured  within  one  fiber  diameter  from  the  surface.  Therefore,  the  change  of  the  a.xial 
stress  and  consequently  the  magnitude  of  the  interfacial  shear  (equation  10)  in  the  real  sys¬ 
tem  should  be  smaller  than  that  calculated.  This  further  suggest  that  the  interface  may  be 
debonded  for  a  distance  about  one  fiber  diameter  from  the  surface.  However,  this  debonding 
does  not  seem  to  change  the  stress  field  significantly,  so  that  the  fiber/matrix  system  still 
behaves  in  an  elastic  manner.  This  is  consistent  with  the  fact  that  the  friction  is  sufficiently 
large  in  this  system,  due  to  asperities  at  the  interface  [4.20),  that  the  stress  distribution  is 
similar  to  a  system  with  perfect  bonding  [3].  Therefore,  we  believe  that  the  residual  stresses 
determined  for  the  sapphire/7-TiAl  system,  especially  in  the  interior  of  the  composite,  are 
not  significantly  affected  by  possible  debonding  near  the  surface. 

6  Conclusions 

By  utilizing  the  optical  fiuorescence  from  chromium  impurities  in  sapphire  fibers,  we  have 
determined  the  residual  stress  distributions  in  embedded  fibers  in  both  a  7-Ti.Al  matri.x  and 
a  polycrystalline  .AI2O3  matrix  composite.  The  a.xial  and  radial  components  of  the  residual 
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stress,  produced  by  differential  thermal  contraction  on  cooling  after  fabrication,  in  the  fibers 
are  determined  as  a  function  of  depth  below  the  composite  surface.  For  both  composite*', 
the  two  stress  components  vary  over  a  distance  of  about  the  fiber  diameter  from  the  surface 
and  then  are  almost  independent  of  distance  further  into  the  composite. 
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TABLE  I 


PROPERTY 

Sapphire 

7-TiAl 

AI2O3 

E  (GPa) 

434 

173 

250 

V 

0.27 

0.33 

0.27 

a  (10-V“C) 

S.3  (a),  9.0  (c) 

11  ~  13 

8.3  <  a  <  9.0 
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Figure  Captions 

Figure  1.  Schematic  diagram  of  the  optical  arrangement  used  to  collect  fluorescence  from  a 
region  in  an  embedded,  transparent  fiber  below  its  surface.  The  stress  induced  shift  in 
frequency  of  a  characteristic  luminescence  line  is  used  to  determine  the  stress  in  the  probed 
region. 

Figure  2.  Illustration  of  the  effect  of  the  finite  depth  of  field  of  the  optical  system  used. 
When  the  excitation  laser  is  focused  at  a  distance  z  below  the  top  surface,  only  the  signal 
generated  in  a  small  volume  abo%e  and  below  the  focal  plane  is  effectively  collected. 

Figure  .3.  The  depth  of  field  function  of  the  microprobe  used.  It  describes  the  relative 
collection  efficiency  as  a  function  of  distance  from  the  focal  plane.  Details  are  given  in  the 
text. 

Figure  4.  The  fluorescence  line  shift  in  the  region  between  the  tension  surface  (z  =  0)  and 
the  compression  surface  (z  =  -513  ^m)  of  a  ruby  beam  under  4-point  bending.  The  solid 
line  representing  the  true  shift  is  obtained  by  deconvoluting  the  measured  shift  with  the 
depth  of  field  function. 

Figure  5.  The  line  shift  depth  profiles  (z  =  0  at  the  surface)  in  a  c-axis  sapphire  fiber 
embedded  in  the  7-TiAl  matrix. 

Figure  6.  The  line  shift  depth  profiles  (z  =  0  at  the  surface)  in  a  c-axis  sapphire  fiber 
embedded  in  a  polycrystalline  .\l2O3  matrix. 

Figure  7.  .\  comparison  of  the  calculated  and  the  measured  normalized  shifts  as  a  function 
of  depth  for  the  sapphire/7-Ti.\l  system. 
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Figure  8.  The  radial  ^Lnd  axial  stresses  as  a  function  of  depth  in  the  sapphire/7-TiAi 
system  determined  from  the  fluorescence  data. 

Figure  9.  .A.  comparison  of  the  calculated  and  the  measured  normalized  shifts  as  a  function 
of  depth  for  the  sapphire/.AIj03  system. 

Figure  10.  The  radial  and  axial  stresses  as  a  function  of  depth  in  the  5apphire/.Al203 
system  determined  from  the  fluorescence  data. 

Figure  11.  The  interfacial  shear  stresses  in  the  sapphire/7-TiAl  and  the  sapphire/ .Ai203 
systems  calculated  from  the  finite  element  computation  and  the  fluorescence  data. 
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Abstract 

The  concentration  dependent  frequency  shift  of  the  fluorescence  from  chromium 
ions  in  sapphire  is  shown  to  be  a  direct  probe  of  the  image  stress  produced  by  the 
chromium  substituting  for  the  smaller  aluminum  ions  in  sapphire.  Our  analysis  is 
based  on  a  comparison  of  the  frequency  shifts  produced  when  the  volume  of  a 
sapphire  crystal  is  changed  by  an  external  pressure  to  that  produced  by  a 
concentration  of  substitutional  solute.  The  volumetric  strain  sensed  by  the 
fluorescing  chromium  ion  when  the  lattice  is  expanded  by  substitutional  solute  is 
approximately  one-third  of  that  determined  from  X-ray  lattice  parameter 
measurements,  in  accord  with  the  Eshelby  analysis  of  the  strain  field  due  to  a 
concentration  of  elastic  defects  in  a  finite  body. 


1  Introduction 


The  question  as  to  the  lattice  distortions  produced  in  a  crystal  by  a  substitutional  solute 
is  of  long-standing  interest,  and  has  been  the  subject  of  considerable  intellectual  inquiry. 
Although  it  has  been  recognized  throughout  that  the  only  proper  description  is  one  that 
takes  account  of  the  atomic  nature  of  the  solute  and  the  appropriate  inter-atomic  potentials, 
early  investigators  did  not  have  available  to  them  this  information.  In  the  absence  of  such 
information,  it  was  necessary  to  address  the  question  using  continuum  descriptions  of  the 
solute  and  its  effect  on  the  lattice.  One  approach,  that  might  be  termed  the  elastic  defect 
model,  was  to  consider  a  solute  as  a  center  of  dilatation  supplemented  by  point  singulaurities 
of  higher  order.  This  enabled  the  methods  of  continuum  elcisticity  pertaining  to  centers  of 
pressure  to  be  employed  in  finding  solutions.  For  example,  Bilby  modeled  an  interstitial  atom 
as  a  misfitting  elastic  sphere  forced  into  a  smaller  hole  in  a  material,  and  derived  solutions 
for  the  interaction  of  a  solute  with  a  dislocation  (Ij.  Later,  Eshelby,  in  one  of  the  classic 
papers  in  materials,  provided  a  rigorous  foundation  for  the  elastic  defect  model  by  deriving 
the  forces  acting  on  an  elastic  singularity  in  a  solid  [2].  In  subsequent  seminal  papers  [3-5], 
Eshelby  set  forth  a  comprehensive  framework  for  the  computation  of  elastic  displacements, 
strains  and  stresses  resulting  from  the  presence  of  point  defects,  inclusions,  precipitates 
and  transforming  regions.  By  envisaging  a  series  of  hypothetical  cutting,  transforming  and 
welding  operations,  and  the  forces  and  displacements  involved,  Eshelby  formulated  analytical 
expressions  for  the  stresses  and  strains  created  by  a  region  transforming  but  constrained  by 
a  surrounding  matrix. 

Two  concepts  introduced  by  Eshelby  are  of  particular  importance.  One  is  that  inside  any 
finite  body  containing  defects  there  exists  an  additional  displacement  field,  and  by  extension 
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strain  and  stress  field,  that  arises  from  the  boundary  condition  that  at  the  surface  of  the 
body  the  normal  stress  must  be  zero.  By  analogy  with  electromagnetic  fields,  Eshelby  termed 
this  additional  elastic  field  the  image  field.  The  image  field  is  independent  of  position  inside 
the  body  and  independent  of  the  size  of  the  body.  The  second  key  concept  was  that  of  the 
“stress-free’’  strain,  the  shape  change  of  the  transforming  region  that  would  have  occurred 
in  the  absence  of  any  constraint.  (This  same  strain  has  been  denoted  as  the  “esgenstrain" 
by  Mura  [6],  It  also  appears  as  the  “natural”  strain  in  recent  literature  on  the  virtual  crystal 
approximation  (VGA)  approach  to  electronic  structure  calculations  [7].)  Using  the  concept 
of  a  stress-free  transformation  strain,  Eshelby  was  able  to  derive  general  expressions  for  the 
strains  produced  in  the  matrix  as  a  result  of  the  transformation  ais  a  function  of  the  shape 
of  the  dlipsoidal  inclusion.  The  matrix  strain  is  shown  to  consist  of  two  components,  the 
strain  n  the  matrix  if  it  were  of  infinite  extent  and  the  image  strain. 

The  existence  of  two  components  to  the  strain  field  in  a  body  containing  defects  naturally 
raises  the  question  a5  to  which  strains  are  actually  being  determined  when  measuring  the 
lattice  parameter  by  diffraction  methods.  .Although  there  was.  for  a  period,  considerable 
confusion  concerning  this  question,  it  now  appears  resolved  that  Bragg  diffraction  methods, 
using  y  -rays  and  neutrons,  provide  a  direct  measure  of  the  spatially  averaged  total  strain 
in  the  body.  (Diffuse  scattering  provides  additional  information  on  the  distribution  of  the 
lattice  distortions  but  does  not  include  information  about  average  lattice  strain.)  Since 
the  image  strain  was  shown  by  Eshelby  to  be  a  constamt  fraction  of  the  total  strain,  the 
magnitude  of  the  image  strain  has  been  deduced  from  the  concentration  dependence  of  the 
lattice  parameter.  However,  there  have  not  been,  as  far  as  the  authors  are  aware,  any  direct 
measurements  of  the  image  strain. 

In  the  course  of  rather  different  studies  [8-10],  directed  towards  the  application  of  optical 
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fluorescence  methocis  to  the  measurement  of  stresses  in  composite  materials,  we  have,  we 
believe,  found  that  the  fluorescence  from  substitutional  solutes  of  chromium  directly  probes 
the  image  stress  in  single  crystals  of  sapphire  (AI2O3).  In  this  short  paper,  we  describe 
how  we  reached  this  conclusion.  It  is  first  necessary,  however,  to  re-examine  the  data  in 
the  literature  on  the  chromium  solid  solution  in  aluminum  oxide,  together  with  our  own 
experimental  data. 

2  Chromium  Substitution  in  Sapphire 

Solid  solutions  of  chromium  in  aluminum  oxide  provide  a  rich  data  base  to  re-examine 
the  question  of  the  short-range  and  long-range  distortion  induced  by  the  substitution  of  a 
mis-fitting  solute.  In  the  literature  there  exist  data  on  the  chromium  concentration  depen¬ 
dence  of  density  [11],  lattice  parameter  determined  by  X-ray  diffraction  [11],  and  frequency 
of  characteristic  optical  fluorescence  lines  [12].  In  addition,  both  lattice  parameter  and  flu¬ 
orescence  line  shifts  as  a  function  of  hydrostatic  pressure  have  been  reported  for  selected 
chromium  ion  concentrations  [13-15].  The  system  .\l2O3-Cr2O3  also  has  the  advantage  that 
there  is  complete  solid  solution  between  the  two  oxide  phases  [16]. 

The  substitutional  nature  of  the  chromium  ion  has  been  confirmed  by  density  and  X-ray 
measurements  (below).  This  is  consistent  with  both  the  existence  of  the  continuous  solid 
solution  with  CrjOa,  the  similarity  in  ionic  radius  between  Cr^"*'  (0.064  nm)  and  (0.057 
nm),  and  the  fact  that  the  corundum  structure  is  a  rather  close-packed  structure  and  does 
not  contain  any  interstices  as  large  as  the  chromium  ion.  Structure  analysis  as  a  function 
of  hydrostatic  pressure  further  shows  that  the  corundum  structure  of  sapphire  compresses 
uniformly,  with  all  the  distances  decreasing  equally  and  all  the  interatomic  angles  remaining 
unchanged  [17]. 
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Chroniium  in  sapphire  is  also  a  chromophore,  giving  rise  to  the  characteristic  red  fluo¬ 
rescence  of  ruby.  The  fluorescence  originates  from  electronic  transitions  in  the  chromium 
ion  produced  by  appropriate  excitation.  Electrons  excited  from  their  ground  state,  '‘Aj,  can 
lose  their  excess  energy  by  a  variety  of  radiative  (and  non-radiative)  processes.  The  most 
pronounced  transitions  being  the  R-fluorescence  lines  at  1.790  and  1.794  eV  in  the  red.  {A 
recent  paper  tabulates  al\  the  radiative  transitions  in  the  visible  part  of  the  spectrum  [18]). 
The  chromium  ion  can  substitute  for  any  of  the  aluminum  ions  in  the  corundum  crystal 
structure  of  sapphire  with  the  result  that  they  are  octahedrally  coordinated  to  neighboring 
oxygen  and  aluminum  ions.  The  octahedron  is  not  regular,  but  elongated  along  the  c-axis 
of  the  lattice,  resulting  in  splitting  of  the  R-fluorescence  lines  to  form  the  doublet  of  Ri  and 
R2  lines.  X-ray  refinements  [19.20],  in  highly  doped  (4a/o  Cr)  sapphire  show  the  chromium 
ion  to  be  slightly  displaced  along  the  c-axis  from  the  position  normally  occupied  by  the  alu¬ 
minum  ion  so  that  the  chromium  site  symmetry  is  further  reduced,  but  nevertheless  retains 
a  three  fold  rotational  symmetry  about  the  c-axis.  The  effect  of  an  applied  stress  is  to  distort 
the  crystal  field  surrounding  the  chromium  ion,  change  its  potential  energy  and  thereby  alter 
the  energies  of  the  radiative  transitions. 

2.1  Optical  Fluorescence  Data 

The  frequencies  of  the  Rj  and  Rj  ,  as  well  as  the  e.xtended  neighbor  fluorescence  N  lines, 
have  all  been  reported  to  shift  with  pressure  and  with  chromium  concentration.  The  pressure 
dependence  is  well  known  [13,14]  and  continues  to  be  subject  to  further  refinement  because 
of  its  central  importance  as  a  pressure  gauge  for  high  pressure  studies.  The  reported  pressure 
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dependence  of  the  frequency  shift  Ai/  may  be  expressed  as: 

Av  = —(21111 -r  n33)P  (1) 

where  fl.j  is  the  piezo-spectroscopic  coefficient.  The  original  work  of  Formatn  gjL.  aL  [13] 
is  reproduced  in  figure  la  for  the  Ri  line  with  the  sum  2nu  +  having  a  value  of  7.7 
cm'^GPa"^  The  results  for  R2  line  is  very  similar  and  therefore  not  included  here.  Recent 
work,  examining  the  fluorescence  from  beams  loaded  in  four  point  bending,  has  shown  that 
the  equation  is  equally  valid  for  tension  [8,9]. 

The  chromium  concentration  dependence  of  the  fluorescence  energy  has  been  the  subject 
of  just  one  previous  study  [12|.  Our  own  results  have  proven  to  be  essentially  identical 
to  those  of  Kaplyanskii  e^:.  ^  [12],  which  are  reproduced  in  figure  lb.  The  data  may  be 
expressed  by  the  linear  dependence: 


Au  =  99c^  (2) 

where  Su  is  measured  in  wavenumber  and  is  the  chromium  concentration  (in  w/o). 

2.2  X-ray  Lattice  Parameter  Measurements 

X-ray  lattice  parameter  measurements  have  been  made  as  a  function  of  pressure  on  rubies 
with  a  number  of  different  chromium  concentrations  [15,17,18].  These  studies  indicate  an 
isotropic  compressibility  with  the  interatomic  distances  decreasing  equally  with  pressure 
and  the  interatomic  angles  remaining  unchanged  up  to  at  least  4.6  GPa  [17].  The  results  of 
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d’ Amour  gLaL  [i5]  is  shown  in  figure  Ic,  which  can  be  fitted  to  the  following  linear  relations: 

a  =  0.47664(1  -  0.001173  P) 

and 

c=  1.3007(1  -0.001332  P)  (3) 

where  the  lattice  pwameters  a  amd  c  axe  measured  in  nanometers,  and  P  is  measured  in 
GPa.  They  also  found  that  the  density  increased  with  increasing  pressure  according  to  the 
Murnaghan  equation  of  state. 

The  X-ray  lattice  parameter  measurements  were  also  made  as  a  function  of  chromium 
concentration  [11]  (figure  Id),  over  a  range  of  chromium  concentration  up  to  2  w/o  by  Jan 
et.  al..  They  indicate  that  the  lattice  parameters  and  derived  density  are  given  by  the  linear 
relationships: 

a  =  0.47591(1  +0.052Tc,n) 
c  =  1.29894(1  +0.0452  c,„) 

and 

px  =  3.9864(1  +0.340  c,„)  (4) 

where  the  density  px  is  measured  in  g/cm^. 

2.3  Density  Data 

Jan  et,.  ^  have  also  made  precision  density  measurements  [11]  using  the  picnometry 
technique.  Their  results  indicate  a  linear  dependence  of  the  density  on  the  chromium  con- 
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cent  rat  ion  given  by  the  relationship: 


p  =  3.9860(1  +0.341c„)  (5) 

The  lineau'  dependence,  as  cleairly  shown  in  figure  2,  is  consistent  with  chromium  being  a 
substitutional  solute  [21].  The  physical  density  determined  above  is,  within  experimental 
error,  the  same  as  the  X-ray  density  measurements  by  Jan  gL.  aL  in  equation  4.  This 
demonstrates  the  equivatlence  of  the  solute  induced  macroscopic  strain  and  the  mean  strain 
measured  by  X-ray. 

2.4  Frequency  Shift  as  a  Function  of  Mean  Strain 

From  the  results  presented  in  figure  1,  the  frequency  shift  cam  be  related  to  the  mean 
volumetric  strain  d^VjV  by  using  tne  following  relationship: 

SV  _  ^  ila  Ac 
‘  V  “  a  ^  c 

Therefore,  from  equations  1  and  3,  amd  equations  2  and  4,  we  obtain  the  following  expressions 
for  the  fluorescence  shifts; 

Ai/^  =  2100  — 

=  660  ^  (6) 

depending  on  whether  the  strain  is  produced  by  an  externaJ  pressure  or  by  solute  doping. 
The  superscripts  P  and  D  denote  pressure  and  solute  doping  respectively. 
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3  Analysis  and  Dbcussion 


The  intriguing  feature  of  the  data  assembled  in  the  previous  section  is  that  the  apparent 
strain  sensed  by  a  fluorescing  chromium  ion  is  different  depending  on  whether  the  lattice 
is  strained  by  the  application  of  a  hydrostatic  pressure  or  internally  by  the  addition  of 
substitutional  solute.  Specifically,  for  the  same  lattice  strain  measured  by  changes  in  the 
lattice  parameter  by  X-ray  diffraction  the  frequency  shift  of  the  chromitim  line  produced  by 
doping  is  about  one-third  (equation  6)  that  produced  by  an  applied  pressure: 

Ax/^  ~  -  Ai/^  (7) 

VVe  believe  that  the  different  pressure  and  solute  concentration  effects  may  be  reconciled 
by  postulating  that  the  frequency  shift  of  the  fluorescence  lines  is  sensitive  to  the  image  strain 
whereas  .X-ray  diffraction  monitors  the  average  lattice  strain.  To  prove  this,  we  proceed  to 
perform  a  set  of  imaginary  experiments  (figure  3),  assuming  for  simplicity  that  the  medium 
is  elastically  isotropic  and  that  the  concentration  of  solute  is  sufficiently  small  that  any  solute 
interaction  effects,  such  as  dipolar  or  higher,  can  be  ignored. 

We  start  with  a  body  of  sapphire  of  volume  V,  free  of  external  pressure  (P  =  0),  containing 
only  a  few  Cr^'*’  ions  to  sense  the  straun,  so  that  the  solute  concentration  /  — ►  0.  The 
frequency  of  the  fluorescence  line  is 


u  =  uiP  =  QJ~*0)  .  (8) 

This  frequency  is  expected  to  be,  in  the  complicated  manner,  dependent  on  the  local  strain 
in  the  vicinity  of  each  Cr^'*'  ion  as  it  is  a  substitutional  solute,  larger  than  the  .A.F'*'  ion  it 
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replaces.  Fortunately,  by  choosing  this  state  as  the  reference  state,  a  knowledge  of  the  effect 
of  this  local  strain  is  unimportant  since  it  is  the  same  for  both  the  state  of  interest  and  the 
reference  state  and  hence  cancels  out. 

In  the  first  experiment,  a  hydrostatic  tensile  stress  —P  is  applied  to  the  body,  so  that 
its  volume  increases  by  AV.  The  strztin  dependence  (due  to  an  applied  stress  —P)  of  the 
frequency  shift  [9]  for  an  isotropic  material  can  be  in  general  expressed  as 

Ai/  =  Ac„  (9) 

where  the  convention  of  summatioii  of  repeated  indices  is  assumed  and  where  A  is  the  piezo- 
spectroscopic  volumetric  strain  coefBcient. Therefore,  the  frequency  shift  due  to  the  stress 
—  P  is  simply 


Ai/^  =  i/(-P,/-.0)  -  i/(P  =  0./-0) 

= 


(10) 


Alternatively,  in  the  second  experiment,  we  can  produce  such  a  volume  change  by  in¬ 
creasing  the  number  of  Cr^  ions  in  the  body  from  virtually  zero  in  our  reference  state  to 
.V,  distributed  at  random.  Following  Eshelby,  the  displacement  field  generated  in  an  infinite 
body  by  ail  the  substitutional  ions  may  be  written  as: 


(r-)  =  c£ 


(11) 


where  C  is  the  dilational  strength  of  an  individual  Cr^"^  ion,  proportional  to  its  misfit. 

In  a  finite  crystal,  the  displacements  are  modified  by  the  existence  of  the  image  term. 
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required  so  that  the  tractions  at  the  external  surface  sue  zero.  The  additional  displacement 
at  position  f  produced  by  the  image  force  is  [3]; 


8x(1-2^)C.V. 

u  (r)  =  ■  — .  . r 

^  ^  3(1 +cr)  V 


12) 


where  <t  is  Poisson’s  ratio.  The  total  displacement  is  then  the  sum  of  these  two  displacements: 


r-fl 


.-tota{/::«\  _  '  '  n  .  8^(1  2(T)  CN  _ 

“  3(1+.) 


(13) 


The  volume  change  produced  by  the  solute  is  then  simply  the  integral  of  the  displacement 
field  over  the  surface  of  the  bodv: 


AV 


Jv 


•  ds 


=  4'C.viHi — 11  =  4,tC.V7 

1  -f  <7 


(14) 


From  this  volume  change,  the  mean  strain  in  the  body  produced  by  solute  doping  is: 


AV  ^  -V 

=  -rr  =  4jr7f7  77  =  4-7C / 


(15) 


where  /  is  the  concentration  of  Cr^'*'. 

However,  the  3tr2un  field  sensed  by  an  individual  substitutional  ion  is  different  from  that 
described  above  as  we  now  show.  Consider  a  smadl  but  finite  volume  element,  W,  of  arbi¬ 
trary  size  but  which  does  not  contain  any  solute  ions  other  thaii  a  single  probe  ion  whose 
fluorescence  we  may  use  to  monitor  the  state  of  this  small  volume  element.  The  volume 
expansion  (due  to  doping  in  the  body)  of  this  element  of  volume  may  be  written  as  before 
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as: 

AW  =  [  (f)  •  ds  (16) 

Jw  Jw 

In  this  equation,  the  first  term  gives  the  contribution  of  the  total  displacement  field,  which 
includes  the  image  term  due  to  all  substitutional  ions  (equation  12),  and  the  displacement 
field  due  to  the  probing  ion  itself,  since  the  displacements  due  to  other  ions  in  equation 
11  not  included  in  W  do  not  contribute  to  the  integration.  Because  the  displacements 
generated  by  the  probing  ion  are  aJready  in  the  reference  state,  to  which  the  volume  change 
is  being  compared,  its  contribution  is  subtracted  by  the  second  term  in  the  equation  16.  So, 
effectively,  only  the  image  displacement  field  due  to  the  presence  of  all  substitutional  ions  in 
the  body  contributes  to  AW.  VVe  may  then  write: 

Air  =  wcf  (17) 

1  -f-  (7 

Therefore,  the  strain  sensed  by  the  solute  ion  is  simply  the  image  strain: 

AW 

^:r  =  ^  =  ^^(7-i)c/  (IS) 


Comparing  with  equation  15,  the  corresponding  frequency  shift  of  the  fluorescence  from  the 
sensing  solute  ion  would  then  be: 


Au^  =  i/(P  =  0,/=^)  -  t'(P  =  0,/-0) 


“V  —  1  7 

A£:r  =  A  =  X  ~ 


1  Ar 

r 


(19) 


If  one  assumes  that  the  Poisson's  ratio  is  1/3,  then  a  comparison  of  equations  10  and  19 
indicates  that  the  effect  of  solute  concentration  on  the  fluorescence  shift  is  about  1/3  that 


11 


produced  by  an  external  pressure,  ass  found  experimentally  (equation  6).  Since  in  the  case 
of  solute  concentration,  the  strain  affecting  the  shift  is  entirely  due  to  the  image  term,  we 
therefore  conclude  that  the  fluorescence  probes  the  image  strain  whereas  X-ray  diffraction 
provides  a  measurement  of  the  total  strain. 
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Figure  Captions 

Figure  1.  Data  on  the  chromium  solid  solution  in  aluminum  oxide,  (a)  Pressure 
dependence  of  frequency  shift  for  the  Ri  fluorescence  line  [13].  (b)  Concentration 
dependence  of  frequency  shift  for  the  Ri  fluorescence  [12]  line,  (c)  Pressure  dependence  of 
the  lattice  parameters  of  sapphire  [loj.  (d)  Concentration  dependence  of  the  lattice 
parameters  of  sapphire  [11]. 

Figure  2.  Linear  dependence  of  the  density  of  alumina  oxide  on  the  chromium 
concentration  [11]. 

Figure  3.  Illustration  of  a  set  of  imaginary  experiments  to  produce  a  volume  change  by 
either  external  stress  or  solute  doping.  The  region  IV  is  small  but  finite  volume  of  arbitrary 
size  that  contains  a  single  probe  ion  used  to  monitor  the  strain. 
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